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ABSTRACT 


Studies  of  uranium-lead  systematic s  in  eight  plutonic  rocks 
indicate  that  the  western  granodiorite  which  intrudes  the  Yellowknife 
volcanic s  has  whole  rock  and  mineral  isochron  ages  of  2700  ±  85  m.y. 
and  2550  ±  70  m.y.  respectively.  Some  new  sulfide  and  galena  data  from 
the  Yellowknife  volcanic s  and  the  sediments  suggest  a  mineralization 
age  of  close  to  263O  m.y..  Fractional  removal  of  lead  from  one  south¬ 
east  granodiorite  rock  sample  yields  a  mineral  isochron  age  of  2660  ± 
45  m.y..  There  is  a  strong  indication  that  the  most  acidic  phase  of 
the  Stock  Lake  intrusive  was  remobilized  at  about  1880  t  75  m.y.. 

The  volatilization  technique  with  an  induction  furnace  was  used 
for  extractions  of  most  of  the  lead  samples  studied  in  this  thesis  and 
this  technique  was.  found  to  be  highly  satisfactory. 

Some  experimental  data  on  mass  discrimination  suggest  that 
different  positioning  of  the  sample  filaments  between  runs  may 
produce  appreciable  mass  dependent  errors  in  the  measured  isotopic 
ratios. 


IV 


ACKNOWIEDGMENTS 


The  writer  wishes  to  express  his  sincere  appreciation  to  Dr.  G. 

L.  Cumming  for  his  enthusiastic  supervision  and  constant  assistance 
which  made  this  thesis  possible.  The  writer  is  also  deeply  indebted  to 
Dr.  H.  Baadsgaard  of  the  Department  of  Geology  who  helped  the  writer 
so  often  and  in  so  many  ways  with  his  skilled  analytical  techniques. 

Lead-free  water  and  other  repurified  acids  for  lead  extraction 
purposes  were  made  available  to  the  writer  by  the  Department  of 
Geology  through  the  courtesy  of  Dr.  Baadsgaard  and  his  associates 
and  is  very  much  appreciated. 

The  automatic  digital  mass  spectrometer  designed  by  Dr.  Cumming 
was  maintained  during  this  project  by  the  able  technicians  Mr.  M. 
Burke,  Mr.  H.  McCullough  and  Mr.  W.  Williams.  Considerable  improve¬ 
ments  in  the  design  of  the  cold  finger  and  vacuum  evacuation  system 
used  for  volatilization  were  made  by  Mr.  D.  Kristie. 

Financial  assistance  was  provided  by  the  University  of  Alberta 
in  the  form  of  a  Graduate  Teaching  Assistantship  (1968-69)  and  a 
University  of  Alberta  Dissertation  Fellowship  (1969-70)  and  is  greatly 
appreciated. 


v 


TABLE  OF  CONTENTS 


CHAPTER  PAGE 

I.  INTRODUCTION . 1 

Introduction  . . 1 

Sample  locations . 2 

Radiometric  dating  .  8 

Mass  spectrometry . 9 

II.  DETERMINATION  OF  U,  Th  AND  Pb . . .  10 

Determination  of  U  and  Th . 10 

Neutron  activation . 11 

Comparison  of  chemical  and  neutron  activation 

methods  for  U  and  Th  determinations  .  13 

Method  used  for  U  and  Th  determinations  .  14 

Extraction  of  Pb  from  rock  sample  .  16 

Extraction  of  Pb  by  volatilization .  17 

Mass  spectrometric  measurement  and 

data  reduction  . .  18 

Least  squares  polynomial  fitting  . 19 

Internal  errors  and  external  reproducibility..  26 
Test  runs  on  standard  samples  and  mass 

dependent  errors  . 28 

III.  VOLATILIZATION  . .  33 

Introduction  .  33 

Furnaces  .  34 


vi 


. 


Analytical  procedures  .  37 

Loading  PbS  on  the  filament  .  38 

Blanks  . . .  39 

IV.  MASS  DISCRIMINATION  IN  A  MASS  SPECTROMETER .  42 

Introduction  . . 42 

Measuring  error  and  fractionation  due  to 

the  positioning  of  the  filament  .  44 

Focusing  effects  on  the  isotope  ratios  .  45 

Conclusions  on  the  mass  discrimination 

experiments  .  50 

V.  WHOLE  ROCK  DATA  AND  U/Pb  AGES  .  53 

Whole  rock  Pb,  U  and  Th  determinations  .  53 

Uranium-lead  radiometric  method  .  56 

Whole  rock  lead-lead  and  uranium-lead  ages  . .  68 

Possible  three  or  higher-stage  system  .  73 

Second  whole  rock  ratios  .  74 

VI.  ISOTOPIC  COMPOSITION  OF  FRACTIONALLY  VOLATILIZED 

LEAD  AND  THE  TEMPERATURE  ISOCHRON . 78 

Introduction  .  78 

Temperature  isochron . 82 

Two-stage  model  temperature  isochron  ages  ...  98 

Mineral  and  whole  rock  isochrons  .  100 

Some  implications  . 105 


Vll 


Temperature  isochron  of  southeast 

granodiorite  .  107 

VII.  MINERAL  LEAD  COMPOSITION  AND  THE  MINERAL 

ISOCHRON . 110 

Introduction  and  mineral  separation . .  110 

Mineral  lead  composition  and  the  mineral 

isochron  . . Ill 

Sulfide  mineral  isochron  of  the 

Yellowknife  volcanics  . 117 

An  indication  of  the  age  of  the 

Yellowknife  metasediments  .  125 

VIII.  SUMMARY  AND  CONCLUSIONS .  126 

Introduction  .  126 

Interpretation  and  comparison  with 

other  results  . 126 

Mass  discrimination .  129 

*************** 

REFERENCES  CITED  .  131 

APPENDIX  I.  SOME  FORMULAS  ON  VECTOR 

DIFFERENTIATION .  136 

APPENDIX  II.  COMPUTER  PROGRAMS  FOR  DATA  REDUCTION 

BY  IBM  360  .  I37 

APPENDIX  III.  DITHIZONE  EXTRACTION  PROCEDURES  .  162 

VITA  .  l64 


viii 


' 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 


LIST  OF  TABLES 


Description  Page 

Isotopic  ratios  of  "Broken  Hill"  galena  31 

Isotopic  ratios  of  "NBS  Equal  Atom  #982"  lead  32 

Blanks  test  runs  on  combined  volatilization 

and  dithizone  purification  of  lead  40 

Ratios  of  a  lead  sample  continously  run 

for  six  hours  46 

Ratios  of  a  lead  sample  (same  sample  as  that 

used  in  Table  4)  47 

Ratios  of  a  lead  sample  (same  sample  as  that 

used  in  Table  4)  48 

Fractionation  test  of  different  focusing  made 

on  different  isotopes  51 

Whole  rock  lead  isotopic  ratios  54 

Pb,  U,  Th  concentrations  and  U-Pb  ages  55 


Comparison  of  U  and  Th  concentration  determinations 


by  isotope  dilution  and  neutron  activation 
methods  57 

Comparison  of  whole  rock  Pb  ratios  between  the 

first  and  second  heatings  75 

Pb  ratios  of  fractionally  volatilized  lead  83 

Slope  ages  of  the  temperature  isochrons  99 


IX 


14  Ratios  of  the  volatilized  lead  of  the 

southeast  granodiorite  108 

15  Mineral  lead  ratios  of  the  western  granodiorite  112 

16  Slope  ages  and  p  values  of  the  mineral  isochrons  118 

17  Pb  ratios  of  sulfide  samples  from  the 

Yellowknife  volcanics  119 

18  Pb  ratios  of  galena  samples  from  the 

Yellowknife  volcanics  121 


x 


II ST  OF  FIGURES 


Figure 

1.  Geological  map  of  the  Yellowknife  area 

2.  Schematic  geological  model  of  the  Yellowknife 

area 

3.  Sample  locations 

4.  Sample  locations 

5.  Neutron-induced  fission  cross-section  as  a 

function  of  energy 

6.  Time  series  of  a  mass  spectrum 

7.  Polynomial  fitting  of  points 

8.  Time  series  of  raw  and  smoothed  data 

9.  Schematic  diagram  of  the  high  temperarure 

brick  furnace 

10.  Schematic  diagram  of  the  sample  assemblage  in 

place  in  the  coil  of  the  induction  furnace 

11.  The  development  of  lead  in  a  two  stage  system 

and  the  lead-lead  diagram 

12.  The  development  of  lead  in  a  three  stage  system 

and  the  lead-lead  diagram  for  lead  that  has 


Page 

3 

4 

5 

6 

12 

20 

27 

29 

35 

36 

61 


evolved  according  to  a  three  stage  f]=constant 
history  63 

13.  Lead-lead  diagram  for  lead  that  has  evolved 

according  to  a  three  stage  f2=constant 

history  65 

14.  Modified  concordia  plot  for  a  two  stage 

evolution  of  lead  67 

15.  Modified  concordia  diagrams  for  lead  that  has 

evolved  according  to  a  three  stage  f]_  and  f2 
constant  history  69 

16.  Whole  rock  lead  isochron  70 

17.  Modified  concordia  plot  of  samples  from 

western  granodiorite  72 

18.  Whole  rock  lead  isochron  (second  heating)  77 

19.  Temperature  isochron  of  YKWG  #1  90 

20.  Temperature  isochron  of  YKWG  #2  91 

21.  Temperature  isochron  of  YKWG  #3  92 

22.  Temperature  isochron  of  YKWG  #4  93 

23.  Temperature  isochron  of  YKWG  #5  94 

24.  Temperature  isochron  of  YKWG  #6  95 

25.  Temperature  isochron  of  YKWG  #7  96 

26.  Temperature  isochron  of  YKWG  #8  97 

27.  All  temperature  isochron  101 


xii 


28.  Effect  of  metamorphism  on  the  change  in  the 

Sr^/Sr^  with  time  in  a  rock  and  its 

constituent  minerals  103 

29.  Variation  in  Sr^/Sr^  vs.  Rb^/Sr^  for  whole 

rock  and  the  effect  on  the  constituent 

minerals  by  a  metamorphic  event  at  time,  t  104 

30.  Lead-lead  diagram  for  whole  rocks  and 

minerals  showing  a  two-event  history  107 

31.  Temperature  isochron  of  the  southeast 

granodiorite  109 

32.  Mineral  isochron  of  YKWG  #1  114 

33  •  Mineral  isochron  of  YKWG  #4  115 

34«  Mineral  isochron  of  YKWG  #5  116 

35»  Mineral  isochron  of  the  Yellowknife  volcanics  122 

36.  Lead-lead  diagram  of  galena  samples  from 

the  Yellowknife  volcanics  124 

37.  Age  measurements  of  the  Yellowknife  area  128 


xiii 


CHAPTER  I 


INTRODUCTION 


1.1.  Introduction 

As  part  of  a  continuing  investigation  of  the  geological  evolution 
of  the  Yellowknife  area,  this  thesis  presents  the  results  and  inter¬ 
pretations  of  a  study  of  lead-uranium  systematic s  on  several  plutonic 
rock  samples  of  this  area. 

All  the  leads  studied  in  this  thesis  are  extracted  from  rocks  and 
minerals  by  volatilization.  It  was  felt  that  the  volatilization 
technique  for  the  removal  of  trace  amounts  of  rock  lead  for  radiometric 
dating  purposes  may  have  some  advantages  over  other  methods,  and  that 
a  closer  study  of  this  technique  might  be  fruitful  in  an  area  well 
investigated  geochronologically. 

Yellowknife  is  located  on  the  north  shore  of  Great  Slave  Lake, 
approximately  700  miles  north  of  Edmonton.  The  area  surrounding 
Yellowknife  has  been  well  studied  by  investigators  using  various 
radiometric  methods  (Keevil  et  al  (1942),  helium  method;  Folinsbee  et 
al  (1956),  K/Ar;  Folinsbee  (1955),  radiation  dosage;  Cumming  et  al 
(1955),  model  lead  ages  on  galena;  Burwash  and  Baadsgaard  (1962), 

K/Ar;  Van  Breeman  (1965),  Rb/Sr;  Robertson  and  Cumming  (1968),  model 
lead  ages  on  galena;  Green  et  al  (1968),  Rb/Sr;  Green  and  Baadsgaard 
(1971),  U/Pb  in  zircons;  Gates  and  Hurley  (1973),  Rb/Sr  on  dykes). 
Besides  checking  the  volatilization  technique  for  whole  rock  radio- 
metric  dating,  it  is  hoped  that  the  addition  of  a  whole  rock  lead- 
uranium  study  will  yield  a  more  complete  picture  of  the  geological 
evolution  of  the  area. 
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1.2.  Sample  locations 


Figure  1  shows  a  map  of  the  Yellowknife  area.  The  regional  geology 
has  been  described  in  detail  by  Green  (1968)  and  is  not  repeated  in 
this  thesis.  A  schematic  geological  model  of  this  area  is  illustrated 
by  Figure  2  with  some  of  the  relevant  age  measurements  on  the  rocks 
indicated.  The  Yellowknife  group  consists  of  acidic  to  basic  volcanics 
overlain  by  a  thick  sequence  of  sedimentary  rocks  and  is  intruded  by 
three  major  plutonic  bodies.  The  southeast  granodiorite  intrudes  the 
sediments  from  the  southeast,  the  Prosperous  Lake  granite  intrudes  the 
sediments  around  Prosperous  Lake  and  the  western  granodiorite  intrudes 
the  volcanics  from  the  west.  One  of  the  problems  posed  at  the  beginning 
of  the  investigation  was  the  determination  of  the  age  of  the  western 
granodiorite  by  means  of  whole  rock  lead-uranium  systematics  in  order  to 
compare  it  with  the  ages  obtained  by  Green  et  al(l968,  Rb/Sr)  and  others. 

Eight  rock  samples  were  collected  in  the  summer  of  1969  from  the 
body  of  the  western  granodiorite.  These  rocks  have  been  dated  by  the 
whole  rock  Rb/Sr  method  by  Green  et  al  ( 1969 )  who  obtained  an  age  of 
26l0  ±  58  m.y.  (  ^Rb=l. 39x10  ^  yr  in  excellent  agreement  with 
data  on  the  volcanics  themselves.  Some  of  the  samples  were 
collected  from  the  same  localities  as  the  material  used  by  Green  et  al 
(1968). 

Twenty-three  sulfide  and  five  galena  samples  collected  by  A.  P. 

Leech  of  the  Department  of  Geology  and  one  whole  rock  sample  from  the 
southeast  granodiorite  obtained  from  the  collection  in  the  Department 
of  Geology  were  also  analysed. 

The  locations  from  which  samples  were  obtained  from  the  western 
granodiorite  are  shown  in  greater  detail  in  Figures  3  and  4  and  are 
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Figure  1.  Geological  map  of  the  Yellowknife  area  (after  Green  and  Baadsgaard,197l) 


4 


ir\ 

vO 


Figure  2.  Schematic  geological  model  of  the  Yellowknife  area 
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Figure  3«  Sample  locations 
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Figure  4*  Sample  locations 


briefly  described  below; 

Sample  YKWG  #1:  Syenite  (acidic  phase). 

Collected  from  the  island  in  Stock  Lake  (see  map)  on  the  south 
side  of  the  highway  to  Yellowknife  airport. 

Location:  114° 23* W,  62°28'n. 

Sample  YKWG  #2:  Syenite. 

Road  cut  on  the  other  side  of  the  highway  opposite  sample 
YKWG  #1. 

Location:  114°22.9'w,  62°28.2’n. 

Sample  YKWG  #3:  Granodiorite. 

9.6  miles  on  road  west  of  Lake  A  and  Lake  B  (see  map). 

Location:  114O30.2’w,  62°27.S'n. 

Sample  YKWG  #4:  Granodiorite. 

3.1  miles  on  road  west  of  Lakes  A  and  B. 

Location:  114035.2'w,  62°27.6’n. 

Sample  YKWG  #5:  Granodiorite. 

1.6  miles  on  road  west  of  Lakes  A  and  B. 

Location:  114°33.o'w,  62°27.4'n. 

Sample  YKWG  #6:  Pegmatite. 

Collected  at  0.9  miles  on  road  east  of  Lakes  A  and  B. 

Location:  114°29.2,W,  62028.3'n. 

Sample  YKWG  #7:  Granodiorite. 

From  0.8  miles  on  road  north  of  the  road  junction  at  Kam  Lake  ( 
see  map). 

Location:  114°24.l’w,  62°26.4'n. 

Sample  YKWG  #8:  Diorite. 


From  the  west  shore  of  Ryan  Lake  (see  map). 

Location:  114°21.8’w,  62°35.7’n. 

1.3 «  Radiometric  dating 

It  is  well  known  that  the  rate  of  decay  of  radioactive  nuclides 
remains  essentially  unchanged  regardless  of  their  geological  enviro- 
ment  .  If  the  amount  and  rate  of  decay  of  the  radioactive  parent 
nuclide  and  the  amount  of  the  radiogenic  daughter  nuclide  may  be 
determined,  the  time  taken  to  accumulate  the  radiogenic  daughter  may 
be  determined.  There  are  many  naturally-occurring  radioactive  nuclides 
each  of  which  is  characterized  by  its  decay  constant  X,  and  the 
daughter  nuclides  it  generates.  Only  a  few  are  suitable  for  use  as 
nuclear  clocks  for  radiometric  dating  (Baadsgaard,  1965).  The  follow¬ 
ing  decay  schemes  are  commonly  used; 


K40  - 

—  Ar40 

Rb87 

—  Sr87 

U238 

Pb206 

U235 

—  Pb207 

Th232 - 

—  Pb208 

Uranium-lead  geochronometry  and  its  application  to  the  problem  at 
hand  are  described  in  Chapter  V. 

An  accurate  determination  of  the  amount  of  the  radioactive 
parent  as  well  as  the  amount  and/or  the  isotopic  composition  of  the 
accumulated  daughter  product  in  the  sample  is  essential.  The  analy¬ 
tical  procedures  used  for  separation  of  U,  Th,  and  Pb  are  described 
in  Chapters  II  and  III. 


1.4.  Mass  spectrometry 


Mass  spectrometry  plays  an  important  role  in  precise  measuring 
of  the  isotope  abundance  ratios  of  an  element  and  the  methods  of 
isotope  dilution  and  solid  source  mass  spectrometry  make  possible 
the  determination  of  very  small  amounts  of  U,  Th,  and  Pb.  A 
description  of  the  mass  spectrometric  techniques  which  include  data 
reduction  and  a  discussion  of  mass  discrimination  are  given  in 
Chapters  II  and  IV. 


CHAPTER  II 


DETERMINATION  OF  U,  Th  AMD  Pb 


2.1.  Determination  of  U  and  Th 

Total  uranium  may  be  determined  by  a  number  of  methods  which 
include  fluorescence,  photometry,  colorimetry,  polarography  and  alpha 
counting  for  small  amounts,  while  classical  gravimetric  and  titrime- 
tric  procedures  may  be  used  for  samples  high  in  uranium.  Thorium  may 
be  determined  by  colorimetry  for  small  amounts  or  by  titrimetric  and 
gravimetric  procedures  for  large  amounts  (Baadsgaard,  1965).  More 
recently  various  neutron  activation  techniques  have  been  applied 
successfully  to  rocks  containing  trace  amounts  of  U  and  Th  (Amiel, 
1962;  Morgan  and  Lovering,  1963,  1968;  Gale,  1967;  Fisher,  1970). 

The  mass  spectrometric  method  of  isotope  dilution  is  used  most 
often  for  the  accurate  determination  of  small  amounts  of  uranium  and 

pqc  oqq 

thorium  contained  in  the  sample.  U  ^  and  Th  ^  are  utilized  as 

235 

isotope' dilutants  or  spikes,  and  U  in  a  sample  is  computed  from 

pqd  pq  r  pqrt 

the  known  ratio  of  U  ;  /U  D:>  and  the  determined  amount  of  U  ;  , 

232 

Th  contained  in  the  sample  can  be  calculated  from  the  amount  of 
1^30  ac^e£  the  measured  ratio  Th^^/Th^^. 

Before  mass  spectrometric  analysis,  U  and  Th  are  separated  from 
the  rock  sample.  The  process  involves  chemical  dissolution  of  the 
rock  powder  (spiked  or  unspiked)  and  the  use  of  ion  exchange  or  other 
methods  to  recover  uranium  and  thorium  from  solution.  Blanks  are 
always  introduced  through  the  chemical  reagents  used  in  the  separation 
and  an  important  problem  is  to  ensure  complete  solution  of  the  U  and 


10 


. 

' 


Th.  Krogh  (1972)  recently  developed  a  high  pressure  decomposition 
technique  for  the  extraction  of  U  and  Pb  from  zircons  with  very 
small  blanks.  Samples  are  decomposed  at  220°C  with  L$!  HF  in  a  self¬ 
sealing  Teflon  capsule  confined  by  a  stainless  steel  jacket  $  uranium  and 
lead  are  isolated  from  the  solutions  on  a  Teflon  anion  exchange 
column  using  Dowex  1  resin.  Because  of  the  relatively  small  amount  of 
chemicals  required,  the  blanks  are  reduced.  This  high  pressure  decom¬ 
position  technique  could  also  be  applied  to  whole  rock  samples. 

2.2.  Neutron  activation 

A  new  method  for  U  and  Th  analysis  has  recently  been 
developed  by  Gale  (1967)  and  others.  This  method  is  based  on  the 
detection  of  delayed  neutrons  emitted  in  the  decay  of  fission  products 
of  uranium  and  thorium.  Details  of  the  analytical  technique  have  been 
described  by  Gale  (1967)  following  some  previous  studies  and  applica¬ 
tions  of  this  technique  by  Amiel  (1962).  The  sample  is  exposed  to  the 
intense  neutron  flux  inside  a  nuclear  reactor,  and  it  is  the  difference 
in  fission  cross  sectional  dependence  on  neutron  energy  which  make 
possible  the  determination  of  U233,  u238  and  Th282. 

It  can  be  seen  from  Figure  5  that  fast  neutrons  induce  fission 
232  23 S 

principally  in  Th  and  U  while  slow  neutrons  produce  fission  in 
U^33.  Since  the  transuranic  elements  such  as  Np23^,  Pu23<^,  etc. 
are  virtually  absent  in  geological  material  (Gale,  1967),  the  detect¬ 
ion  of  delayed  neutrons  following  exposure  of  a  sample  to  either  slow 
(thermal)  or  fast  neutrons  is  characteristic  of  the  presence  of  U233 
or  Th232  and/or  U238. 


i  -  u  ’  *-  1 J 


0.01  0.05  o.l  0.5  l.o  2  5  10 

Neutron  energy  (Mev) 

Figure  5.  Neutron-induced  fission  cross-section  as  a  function  of 
energy,  (after  Gale,  1967) 
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Two  measurement s,  one  in  a  mixed  flux  which  yields  primarily 

235 

U  fission  products  and  a  second  with  the  slow  neutrons  screened 

out  with  a  Cd  shield  are  required.  The  interference  of  on  Th^^ 

of  the  second  measurement  can  be  corrected  from  the  known  ratio  of 

natural  uranium  and  the  first  measurement.  Gale  (1967)  mentioned  a 

28 

possible  interference  from  the  large  T-ray  flux  emanating  from  A1 
formed  in  the  irradiated  sample  and  noted  that  the  effect  is  especial¬ 
ly  serious  in  feldspar  rich  samples  with  high  aluminium  content. 
Cumming  (personal  communication)  has  shown  that  this  effect  may  be 

controlled  and  that  the  only  serious  interference  is  from  the  reaction 

17  17  17  1 7 

0  '(n,p)N  ,  N  decaying  to  an  excited  state  of  0  '  which  then  decays 

by  neutron  emission,  and  that  this  interference  is  only  important  at 

very  low  U  and  Th  concentrations. 

The  neutron  counts  of  a  sample  are  compared  to  those  of  standards 

and  the  amount  of  U  and  Th  contained  in  the  sample  can  be  calculated. 

During  the  course  of  this  work  it  was  possible  to  compare  results 

from  the  neutron  activation-delayed  neutron  technique  with  isotope 

dilution  measurements. 


2.3 .  Comparison  of  chemical  and  neutron  activation  methods  for 

U  and  Th  determinations 

When  comparing  methods  some  pertinent  questions  to  be  asked  are; 
(i).  Can  the  particular  method  be  conveniently  carried  out  in  one’s 
laboratory?  (ii).  How  much  blank  error  is  involved?  (iii).  How  much 
time  is  required?.  Based  on  these  questions  the  advantage  of  one 
method  over  the  other  may  be  evaluated.  The  chemical  method  is  more 


■ 
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readily  available  than  the  neutron  activation  technique  since  the 
latter  requires  the  use  of  a  nuclear  reactor.  As  there  is  no  need  for 
treatment  of  the  sample,  no  blank  is  involved  in  neutron  activation 
while  blanks  are  always  introduced  in  the  chemical  method,  but  with 
the  use  of  super-pure  chemicals  and  improved  decomposition  techniques 
the  blank  can  be  reduced,  controlled  and/or  corrected.  The  time 
required  for  a  typical  single  analysis  is  quite  long  for  the  chemical 
method  compared  to  that  required  by  the  neutron  activation  technique. 

Perhaps  the  biggest  advantage  of  the  neutron  activation  technique 
over  the  chemical  method  is  that  the  activation  technique  is  non¬ 
destructive.  The  sample  that  had  been  analysed  for  uranium  and  thorium 
can  also  be  analysed  for  lead  thus  eliminating  the  problem  of  sample 
homogeneity.  This  is  very  important  because  large  inhomogeneities  in 
the  distributions  of  lead  and  uranium  have  been  demonstrated  to  have 
occurred  in  meteorites  by  Marshall  (1962)  and  Goles  et  al  (quoted 
from  Gale,  1967)  and  could  also  occur  in  terrestrial  materials. 

2.4.  Method  used  for  U  and  Th  determinations 

Chemical  procedures  followed  in  the  separations  of  U  and  Th  are 
shown  by  the  following  flow  diagram; 


Whole  rock  powder  (~200  mesh), 

about  one  gram. 

. 

Add  U2-^  and  Th2^  spikes  and  dry  at  250°  C 

l  r 


. 
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The  neutron  activation  technique  was  also  used  (  G.  L.  Cumming, 
personal  communication)  for  U  and  Th  determination  on  these  samples, 

A  comparison  of  the  data  from  the  two  methods  is  given  in  Chapter  V. 

2.5.  Extraction  of  Pb  from  rock  sample 

Chemical  extraction  of  rock  lead  involves  a  complete  dissolution 
of  the  rock  powder  and  subsequent  recovery  of  Pb  from  solutions  by 
anion  exchange,  dithizone  extraction  or  other  methods.  Large  amounts 
of  pure  reagents  have  to  be  prepared.  Chemical  separation  of  suffi¬ 
cient  lead  for  analysis  from  rocks  which  contain  trace  amounts  is 
difficult.  Contamination  from  chemical  reagents  and  from  air  borne 
Pb  at  those  low  concentration  levels  may  constitute  an  appreciable 
part  of  the  final  sample  if  care  is  not  taken. 

Air  borne  contamination  may  be  eliminated  to  a  large  extent  in 
a  particle-free  laboratory  in  which  positive  air  pressure  with  respect 
to  its  surrounding  is  maintained  using  filtered  air.  A  similar  local 
particle-free  enviroment  may  also  be  obtained  by  using  laminar  flow 
hoods. 

Purification  of  chemicals  from  commercial  reagent-grade  material 
is  customarily  done  by  redistillation.  Tilton  et  al  (1955)  prepared 
pure  HF  by  passing  hydrogen  fluoride  gas  through  a  filter  of  fine 
Telflon  shavings  and  bubbling  the  gas  into  pure  water  in  an  ice- 
cooled  platinum  vessel.  The  Pb  contamination  was  reported  to  be 
0.0002^g/ml.  of  HF.  Tatsumoto  (1969)  refined  Tilton  et  al’s  (1955) 
method  by  freezing  the  Telflon  filtered  hydrogen  fluoride  gas  in  a 
Kel-F  tank  with  liquid  nitrogen  and  then  redistilling,  refreezing 
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and  subsequently  introducing  the  gas  into  pure  water  in  a  Telflon 
bottle.  Pb  contamination  was  reported  to  be  0. 00008 ^ag/ml.  of  HF. 
Techniques  for  the  preparation  of  highly  pure  chemical  reagents  are 
being  constantly  sought;  Mattinson  (1971  and  1972)  described  a  simple 
two-bottle  subboiling  still  technique  for  purifications  of  HF,  HC1  and 
HNO^.  Measured  lead  contamination  levels  for  HF  and  HC1  were  reported 
to  be  0.002  ppb  to  0.005  ppb,  and  for  HNO^ ,  0,02  ppb. 

The  use  of  highly  pure  chemicals  combined  with  the  improved 
decomposition  technique,  such  as  the  one  described  by  Krogh  (1972), 
the  particle-free  working  enviroment  and  the  high  precision  of  mass 
spectrometric  analyses  will  no  doubt  help  to  improve  the  precision 
and  reproducibility  of  the  isotope  dating  methods. 

For  quantitative  determination  of  Pb,  isotope  dilution  is  often 
used.  A  100  per  cent  dissolution  of  Pb  into  the  solutions  from  rock 
powder  is  customarily  assumed. 


2.6.  Extraction  of  Pb  by  volatilization 

The  volatilization  technique  offers  distinct  advantages  over 
chemical  methods;  air  borne  lead  contamination  during  chemical  pro¬ 
cessing  is  reduced  to  a  minimum  and  it  is  not  necessary  to  prepare 
large  amounts  of  pure  reagents.  Further,  as  shown  by  Cumming  et  al 

(1970)  fractional  removal  of  lead  from  a  rock  sample  by  volatilization 

?D7  on; 

yields  lead  of  varying  isotopic  composition  and  the  ratios  Pb  '/Pb  4 
and  Pb^/Pb204  are  linearly  related.  Its  analogue  to  a  mineral 


isochron  is  demonstrated  in  this  thesis  (see  Chapter  VII )  and  thus  a 
mineral  isochron  may  be  obtained  in  a  simple  fashion  utilizing  the 
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inhomogeneous  distribution  of  volatilised  lead  in  the  rocks  without 
having  to  actually  separate  the  mineral  phases  (which  is  customarily 
done  in  the  chemical  method).  (Perhaps  a  '’mineral"  isochron  of  a  rock 
can  also  be  obtained  by  fractional  removal  of  lead  from  a  rock  sample 
by  chemical  leaching). 

For  quantitative  analyses  of  lead,  the  isotope  dilution  method 
is  used  and  total  evaporation  of  lead  from  the  rock  powder  is  neces¬ 
sary.  It  has  been  demonstrated  by  Starik  et  al  (1957)  and  Marshall  and 
Hess  (i960)  that  the  isotopic  composition  of  the  volatilised  lead  is 
practically  the  same  as  that  of  lead  isolated  chemically.  Black  et 
al  (1971)  estimated  the  volatilization  yield  (at  1250°C  for  4  hours 
under  vacuum)  on  one  U.S.G.S.  standard  sample  previously  analysed 
chemically  by  Tatsumoto  and  obtained  104%  recovery  for  the  volatili¬ 
zation  process  with  respect  to  the  chemical  method.  It  thus  becomes 
apparent  that  total  lead  evaporation  from  rock  is  obtainable  by 
prolonged  heating  at  high  temperature  under  vacuum. 

Purification  of  volatilised  leads  can  be  easily  done  by  the 
dithizone  extraction  technique.  All  the  leads  studied  in  this  thesis 
except  those  from  galenas  were  extracted  by  volatilization.  Details 
of  the  extraction  procedures  and  further  discussions  of  this  technique 
are  given  in  Chapter  III. 

2.7*  Mass  spectrometric  measurement  and  data  reduction 

Since  the  major  analytical  technique  is  mass  spectrometry  and 
the  final  data  depends  on  the  quality  of  the  mass  spectrometer  results, 
it  is  very  desirable  that  the  mass  spectrometric  technique  be  developed 
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to  the  point  where  precise  reproducible  data  can  be  obtained  on  a 
routine  basis.  Cumming  et  al  (1971)  described  an  automatic  digital 
mass  spectrometer  developed  in  this  Department  and  used  by  the  author. 

A  brief  version  of  a  least  squares  polynomial  fitting  technique  for 
the  data  reduction  was  also  given  by  Cumming  et  al  (1971)  and  is 
described  in  greater  detail  in  this  section. 

2.7.1.  Least  squares  polynomial  fitting 

In  determining  the  isotopic  composition  of  an  element  using  a 

t 

mass  spectrometer  it  is  customary  to  make  repeated  scans  over  the  mass 

range  of  interest  and  measure  the  relative  abundances  of  the  isotopes 

often  enough  to  obtain  results  within  useful  statistical  limits.  Errors 

are  easily  made  because  of  long-term  and/or  short-term  drift  due  to 

changes  in  resistances,  temperature  variations  and  even  personnel 

differences  in  measuring  peak  heights  (Kanasewich,  1968).  The  least 

squares  polynomial  fitting  technique  is  aimed  to  minimize  these  errors. 

We  make  M  repeated  measurements  of  a  mass  spectrum  where  there 

are  N  peaks.  An  example  of  three  repeated  scans  of  a  mass  spectrum  of 

four  peaks  with  slightly  growing  signal  is  shown  in  Figure  6a. 

Let  P„,  where  i  =  1,2,3,.... 

j  =  1,2,3,....  be  the  height  of  the  jth 

peak  of  the  ith  scan  in  the  spectrum.  Let  X-»  where  j  =  1,2,3,4»«.« 

3 

be  the  ratio  of  the  first  peak  to  the  jth  peak  in  the  scan.  We  assume 

X-  to  be  constant  throughout  the  measurements. 

3 

Let  the  polynomial  representation  of  the  peak  as  a  function  of 


time  be 


Amplitude  Amplitude 
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Figure  6a.  Time  series  of  a  mass  spectrum. 
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Figure  6b.  Time  series  of  a  normalized  mass  spectrum. 
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P(t)  bQ  +  b^t  +  b^t  +  b^t^  +  ...  +  b^t^  or  simply  as 

Pit)  =  ‘^0bkt  »  where  t  is  the  time  when  the  peak  occurs  in 

the  spectrum. 


Assume  that  the  isotopic  ratios  are  independent  of  time  and  let  R  be 

3 

^ j  ^1  —  ratio  Ox  the  jth  peak  to  the  first  in  a  scan 

and  normalize  the  measured  peak  heights  to  the  height  of  the  first 
peak  in  each  scan  by  the  parameter  x,.  The  time  series  of  the  nor- 

J 

malized  spectrum  will  then  look  like  that  shown  in  Figure  6b. 

The  normalized  peak  heights  are  not  necessarily  of  equal  height 
due  to  the  growth  or  decay  of  signal  during  the  measurements.  It  is 


our  aim  to  solve  for  the  normalization  parameters  x • ,  j 


» 


•  •  •  • 


N  which  will  normalize  the  peak  heights  to  fit  a  single  polynomial. 


Relative  isotopic  abundances  of  the  isotopes  in  a  mass  spectrum  with 
respect  to  the  first  peak  in  this  mass  spectrum  can  then  be  obtained 
by  taking  the  reciprocals  of  the  normalization  parameters. 

Letting  t^  be  the  time  where  occurs  in  the  spectrum  then 
we  can  express  the  measured  and  normalized  peak  heights  as; 
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Regrouping  the  above  equations,  we  have 
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Writing  in  explicit  matrix  form, 
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or  simply  as 


Mb  =  p 

where  M  is  a  (m*n)  by  (n+k)  matrix  with  (m*n)  >  (n+k) 
b  is  a  (n+k)  column  vector 
p  is  a  (m«n)  column  vector. 

The  system  Mb=p  contains  more  equations  than  unknowns  and  in  general 
no  solution  for  b  is  available.  However,  we  are  interested  in  finding 
a  vector  b  that  is  least  in  error.  More  precisely,  if  we  take  a 
vector  b  and  calculate  Mb=q,  say,  we  want  the  discrepancy  between 
q  and  p  to  be  as  small  as  possible.  The  best  overall  criterion  is  to 
minimize  the  sum  of  the  squares  of  the  discrepancies,  which  gives  a 
least  squares  solution  for  Mb=p. 

The  discrepancies  are  given  in  a  vector  (q  -  p),  and  the  sum  of 
the  squared  elements  then  is  given  by  the  square  of  its  norm,  or 
expressed  as 

f=(q-p)T(q-p). 

Using  q=Mb,  we  can  write 

f=  (Mb-p  )T  (Mb-p  ) 

mm  m  m  m  m 

=b  M  Mb-p  Mb-b  M  p+p  p 

rp  rn  rn  m  rn 

=b  M  Mb-2b  M  p+p  p 
T  T  T 

where  p  Mb=b  M  p,  since  both  are  simple  numbers. 

The  problem  then  can  be  redefined  by  saying  that  we  want  to  find  a 
vector  b  for  which  f  has  a  minimum.  Setting  partial  derivatives  to 
zero  (see  Appendix  I), 
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=  2MTMb  -  2MTp  =  0 
Therefore,  M^Mb  =  M^p 

so  that  b  =  (MTM)-1MTp. 

T  —1 

(MM)  is  a  symmetric  matrix,  the  determinant  of  which  should  be 
different  from  zero. 


T  —1  9  p 

The  variance-covariance  matrix  of  b  is  (M  M)  <r  ,  where  <r  is 


estimated  from 


-  \2 


2  IN,  -  v 
O'  =  7?  (p.  -  p.  ) 

N  1=1 

N=nm-(n+k)is  the  number  of  degree  of  freedom,  and  p^  is  the 
normalized  smoothed  value  of  p^.  Thus  the  errors  of  the  ratios  and 
correlation  coefficients  between  ratios  may  be  obtained  directly  from 
the  variance-covariance  matrix.  The  solutions  of  those  equations  are 
obtained  using  programs  from  the  IBM  scientific  subroutine  package. 


2.7*2.  Internal  errors  and  external  reproducibility 

The  size  of  the  errors  is  controlled  by  the  order  of  polynomial; 
although  it  is  possible  to  use  a  high  order  polynomial  equation  whose 
graph  will  pass  each  of  the  data  points  as  closely  as  desired,  it  is 
suggested  that  the  polynomials  be  limited  to  no  more  than  five  terms 
(depends  on  number  of  scans  and  can  be  varied).  The  search  for  an 
exact  equation  presupposes  that  each  of  the  data  points  is  in  itself 
exact,  having  no  error  associated  with  it.  Since  it  is  relatively 
difficult  to  deduce  the  shape  of  a  polynomial  equation  merely  from 
its  coefficients,  the  development  of  such  an  exact  equation  can  pro¬ 
duce  some  rather  implausible  results.  Figure  7  exemplifies  such 


situation. 
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Figure  7.  Polynomial  fitting  of  points.  (Dashed  line:  Lower 
order  of  polynomial;  Solid  line:  Higher  order) 


The  solid  line  represents  a  higher  order  of  polynomial  made  to  fit 
through  all  the  points.  The  path  is  as  questionable  as  the  use  of 
linear  segments  to  connect  successive  values  of  the  points  on  the 
graph,  especially  if  the  set  of  data  contains  one  or  more  points  which 
are  grossly  in  error.  A  more  plausible  alternative  is  the  curve  shown 
as  a  dashed  line  which  has  a  lower  degree  of  polynomial  than  the  solid 
line.  The  internal  errors  are  controlled  by  the  order  of  the  polynomial 
but  the  external  errors  would  probably  be  more  realistic  with  the 
order  of  the  polynomial  limited  to  a  relatively  small  value.  In 
practice  we  make  10  to  15  scans  over  the  spectrum  yielding  40  to  60 
peaks  and  since  there  are  three  unknown  terms  \  even  a  relatively 
high  order  of  polynomial  (10  to  15  degree)  may  be  used  and  still  have 


a  large  number  of  degrees  of  freedom  in  the  least  squares  solution. 


.  •  J  > 

tw  r  o  ?-■  d  ■  ‘ 


2.7.3.  Data  recording  and  reduction 

All  lead  samples  were  run  as  PbS  on  a  single  tantalum  filament  and 
were  recorded  on  magnetic  tape  as  described  by  Cumming  et  al  (1971 ) 
and  processed  by  IBM36O  computer  as  follows; 

(1) .  Each  single  sweep  over  the  mass  spectrum  is  smoothed  by  a  box¬ 
car  filter  with  a  running  average  of  a  number  of  points  specified  by 
the  input  to  reduce  high  frequency  noise. 

(2) .  To  locate  the  times  where  the  peak  heights  occur  in  the  spectrum 
the  data  is  heavily  smoothed  with  a  25  point  running  average  and  the 

times  tn,  t12,  ...  tln,  t21,  t22,...  t2n,  . tmn  are  picked  from 

where  the  maxima  occur  (the  centers  of  the  smoothed  peaks). 

(3) .  The  peak  heights  Pu,  ?u,  ....  P^,  P21>  P22 . Pmn  are 

then  obtained  by  going  back  to  the  smoothed  data  of  (l)  and  take  3  to 
5  points  average  around  each  of  the  maxima.  The  average  of  3  to  5 
points  around  each  of  the  maxima  is  equivalent  to  a  second  smoothing 
with  a  box-car  filter  of  different  length  than  in  (l). 

(4) .  Measured  peak  heights  were  corrected  for  baseline  and  used  for 
polynomial  fitting. 

The  recorded  raw  data  and  the  once  smoothed,  twice  smoothed  and 
heavily  smoothed  data  of  a  scan  of  mass  spectrum  of  a  lead  sample  are 
shown  in  Figure  8. 

The  computer  programs  are  listed  in  Appendix  II. 

2.8.  Test  runs  on  standard  samples  and  mass  dependent  errors 

To  test  the  mass  spectrometric  technique,  the  author  prepared 
10  independent  "Broken  Hill  #1"  and  9  "NBS  Equal  Atom  #982"  standard 
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Figure  8.  Time  series  of  the  raw,  once 
smoothed,  twice  smoothed  and 
heavily  smoothed  data  of  a 
mass  spectrum  of  lead. 
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samples.  Data  were  recorded  on  tape  and  reduced  by  the  polynomial 
fitting  technique  described  above.  Results  are  shown  in  Table  1  and 
Table  2  with  the  error  shown  at  one  standard  deviation. 

Examination  of  the  results  show  external  errors  of  0.1  4$, 

0.21$  and  0.28$  for  ratios  Pb2°6/Pb20/+,  Pb207/Pb20/f  and  Pb208/pb20^ 
respectively  which  represents  0.07$  error  for  every  unit-mass  diffe¬ 
rence  in  the  ratio. 

This  mass  dependent  error,  known  as  mass  discrimination,  is  conce¬ 
ded  by  most  investigators  to  be  a  major  factor  which  might  introduce 
errors  in  the  final  results  (Russell  and  Farquhar,  I960).  It  has  been 
demonstrated  by  Cumming  et  al  (1971),  Cooper  and  Richards  (1966), 

Doe  et  al  (1965)  and  many  others  that  if  the  mass  discrimination 
errors  were  eliminated,  a  very  precise  mass  spectrometric  measure¬ 
ment  (about  0.05$  at  one  standard  deviation  or  better)  could  be 
obtained.  In  this  work  an  empirical  correction  factor  is  applied  to 
all  measurements,  based  on  a  smoothed  mass  discrimination  factor 
obtained  from  the  measured  values  on  the  NBS  Equal  Atom  #982  standard, 
so  the  major  mass  discrimination  effects  should  be  removed  from  the 
data  (Cumming  et  al,  1971). 

In  the  early  stages  of  this  investigation  while  the  author  was 
developing  his  mass  spectrometric  skill,  he  was  bothered  by  this  mass 
dependent  error  and  decided  to  pursue  this  problem  a  little  further. 
Some  studies  of  this  phenomenon  are  given  in  Chapter  IV. 
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Table  1-  Isotopic  ratios  of  "Broken  Hill  #1”  galena 


Average 


Pb206/Pb20/4- 

Pb207/Pb20i+ 

pb208/pb2°4 

16.028 

15.438 

35.805 

.015  * 

.018 

.038 

16.002 

15.403 

35.702 

.028 

.017 

.054 

16.063 

15.471 

35.902 

.  040 

.042 

.110 

16.056 

15.474 

35.928 

.030 

.022 

.044 

16.052 

15.466 

35.902 

.010 

.008 

.019 

16.046 

15.473 

35.941 

.021 

.022 

.045 

16.010 

15.389 

35.630 

.021 

.020 

.060 

16.030 

15.428 

35.786 

.015 

.015 

.035 

16.068 

15.486 

35.935 

.028 

.026 

.049 

16.019 

15.411 

35.719 

.013 

.013 

.030 

16.037  ±  .022  15.444  ±  .033 

(.14/o)  (.21$ 


*  Internal  estimate  of  standard  deviation 


35.825  ±  .107 
( .300) 
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Table  2.  Isotopic  ratios  of  "NBS  Equal  Atom  #982”  lead. 
pb2°6/pb2°4  pb2°7/pb204  Pb208/Pb204 


36.896 

17.283 

37.074 

.072* 

.026 

.037 

36.756 

17.195 

36.838 

.063 

.034 

.056 

36.762 

17.183 

36.842 

.055 

.018 

.047 

36.807 

17.221 

36.946 

.047 

.025 

.068 

36.911 

17.254 

36.978 

.100 

.048 

.106 

36.789 

17.179 

36.793 

.037 

.018 

.039 

36.819 

17.200 

36.864 

.041 

.020 

.043 

36.780 

17.179 

36.836 

.029 

.014 

.031 

36.760 

17.163 

36.754 

.034 

.016 

.036 

Average:  36.809  ±  .054  17.206  ±  .037  36.881  ± 

(.14$)  (. 21 $)  (. 26 $) 


*  Internal  estimate  of  standard  deviation. 
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CHAPTER  III 


VOLATILIZATION 


3 » 1 «  Introduction 

In  this  Chapter  the  volatilization  technique  used  for  lead  extrac¬ 
tion  and  the  loading  of  the  sample  onto  the  filament  are  described. 

The  volatilization  method  has  been  used  for  trace  element  separation 
for  some  time.  Edwards  and  Urey  (1955)  separated  the  alkalies  from 
meteorites  by  volatilization,  Edwards  and  Hess  (1956)  carried  out  pre¬ 
liminary  experiments  for  the  same  kind  of  separation  of  lead  from 
meteorites  and  Starik  et  al  (1957)  investigated  in  great  detail 
different  modes  of  lead  occurrence  in  pitchblende,  monazite,  uraninite 
and  a  series  of  granites  by  volatilization.  Marshall  and  Hess  (i960) 
separated  lead  from  meteorites  by  volatilization  and  dithizone  puri¬ 
fication  with  small  contamination.  Many  other  investigators  including 
Masuda  (1964),  Hamilton  (1965),  Sinha  (1969),  Welke  et  al  (1968), 
Moorbath  et  al  (1968)  and  Cumming  et  al  (1970)  have  employed  this 
technique  for  lead  separation. 

Recently,  the  volatilization  technique  was  used  by  Silver  (1970, 
1971),  Huey  et  al  (1971)  and  Doe  and  Tatsumoto  (1972)  for  separation 
of  lead  from  lunar  samples. 

The  volatilization  procedures  used  in  this  thesis  for  removal  of 
trace  lead  from  rock  are  similar  in  principle  to  that  of  Marshall  and 
Hess  (i960).  Powdered  rock  is  volatilized  in  a  graphite  crucible  inside 
a  furnace  under  vacuum.  Lead,  condensed  on  a  water-cooled  quartz  tube 
(called  a  cold  finger),  is  purified  by  the  dithizone  extraction  method. 
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3.2.  Furnaces 


(A) *  High  temperature  brick  furnace 

At  the  beginning  of  this  investigation  a  high  temperature  brick 
furnace  was  built.  Figure  9  shows  a  sketch  of  the  furnace  and  sample 
assemblage. 

The  furnace  was  heated  by  four  silicon  carbide  elements  which 
surround  the  outer  ceramic  liner.  A  heavy  duty  transformer  was  used 
to  control  the  input  current  to  the  temperature  regulator  which  has 
a  quoted  accuracy  of  +6°C.  Due  to  the  large  temperature  gradient 
along  the  furnace,  the  sample  should  be  placed  in  the  middle  of  the 
furnace  where  maximum  temperatures  occur.  Overall  performance  of  this 
furnace  was  unsatisfactory. 

(B).  Induction  furnace 

After  some  of  the  preliminary  work  was  done  using  the  high 
temperature  brick  furnace,  an  induction  furnace  became  available  and 
produced  highly  satisfactory  results.  Figure  10  shows  a  schematic 
diagram  of  the  sample  assemblage  in  place  in  the  coil  of  the  induct¬ 
ion  furnace.  The  induction  coil  induces  a  large  current  in  the 
graphite  crucible  containing  the  rock  sample  and  heats  up  the  sample 
to  the  desired  temperature  without  excessive  heating  of  the  outer 
quartz  tubing. 

The  temperature  of  the  crucible  is  calibrated  against  the  current 
settings  with  an  optical  pyrometer.  Although  the  exact  temperature 
setting  is  not  important  in  most  of  this  work  (see  Chapter  VI  on 
temperature  isochrons),  the  actual  temperature  is  believed  to  be  close 


■ 
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cold  water  out 


Figure  10.  Schematic  diagram  of  the  sample  assemblage 

in  place  in  the  coil  of  the  induction  furnace. 


to  that  indicated  by  the  pyrometer  (±50°K).  With  an  induction  furnace 
the  desired  temperature  can  be  reached  in  a  few  minutes. 

With  its  easy  operation  and  fast  heating,  the  induction  furnace 
is  much  to  be  prefered  for  this  work  and  most  of  the  author's  lead 
samples  were  obtained  using  this  furnace. 

3«3»  Analytical  procedures 

The  graphite  crucibles,  cold  fingers,  and  the  inner  quartz  tubes 
were  cleaned  by  soaking  them  in  concentrated  hot  nitric  acid  for 
several  hours  and  then  rinsed  with  pure  water  (several  times  redis¬ 
tilled).  The  crucibles  were  then  outgassed  in  the  furnace  to  about 
1400° C  for  approximately  one  hour  before  being  used. 

The  rock  sample  was  first  broken  into  small  fragments  and 
weathered  material  was  discarded.  The  remaining  fragments  were 
crushed  to  100-200  mesh,  mixed  with  a  small  amount  of  powdered 
graphite  (not  essential)  and  placed  in  a  crucible  inside  the  inner 
quartz  tube.  With  the  cold  finger  in  place  slightly  above  the  crucible, 
the  system  was  then  evacuated  and  current  slowly  increased.  For  a 
powdered  sample  the  evacuation  and  the  preliminary  heating  must  be 
very  gentle  otherwise  the  powdered  rock  sample  will  pop  up  and  dirty 
the  walls  of  the  assemblage.  For  a  chipped  rock  sample  the  above 
precaution  is  not  necessary.  After  the  system  has  been  evacuated  to 
10  ^  torr,  the  current  can  be  more  rapidly  increased  to  obtain  the 
desired  temperature. 

Lead  and  other  impurities  deposited  on  the  cold  finger  are 
dissolved  in  5  c.c.  of  redistilled  3N  HNO^  and  purified  by  the 


' 

-  Hi' 


dithizone  extraction  technique  (see  Appendix  III).  All  chemical  ware 
was  initially  cleaned  with  distilled  water  and  placed  in  hot  nitric 
acid  for  several  hours  or  overnight  before  a  final  rinsing  with  pure 
water  and  then  covered  with  parafilm  ready  for  use. 

Purification  of  lead  has  a  very  important  effect  on  the  emission 
of  ions  from  the  lead  sample.  For  best  emission  results,  chemical 
solutions  used  for  purification  of  lead  should  be  freshly  prepared 
and  the  extracted  PbS  be  run  as  soon  as  possible. 

3»4»  Loading  PbS  on  the  filament 

The  loading  of  the  sample  on  the  filament  is  a  critical  step  in 
the  analytical  procedure  and  the  process  used  is  described  below. 

The  tantalum  filament  (carefully  made  to  ensure  that  it  is 
tightly  and  flatly  welded  on  the  filament  posts)  is  outgassed  at 
3*5  amp.  for  3  hours  or  at  3»0  amp.  overnight  under  vacuum  (10-^  torr 
or  better).  The  sample  is  loaded  on  the  filament  in  the  following 
way; 

(1) .  With  a  capillary  tube  (works  best  if  the  narrow  end  is  broken 
off  flat)  pick  up  few  drops  of  liquid  containing  PbS  from  the  centri¬ 
fuge  tube  and  wet  the  center  of  the  filament  with  a  small  drop  of 
liquid.  Excess  liquid  can  be  easily  drawn  up  into  a  spare  capillary 
tube  by  touching  the  liquid  with  the  end  of  the  tube. 

(2) .  Pick  up  some  PbS  and  force  it  down  onto  the  filament  by  touching 
the  wetted  portion  of  the  filament  with  the  end  of  the  capillary  tube. 
The  sample  PbS  will  run  down  onto  the  filament.  Draw  out  excess  liquid 
with  the  spare  tube. 
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(3) .  Apply  a  current  of  up  to  1.5  amp.  through  the  filament  to 
evaporate  the  liquid  drop  containing  PbS.  Decrease  the  current  when 
the  drop  is  about  to  disappear.  Load  more  sample  onto  the  filament  by 
repeating  steps  (2)  and  (3)  until  desired  amount  is  loaded. 

(4) .  Apply  1.5  amp.  of  current  through  the  filament  and  wait  until  the 
sample  is  dried,  then  gently  increase  the  current  to  the  point  at 
which  the  ammonia  fumes  slowly.  It  is  important  not  to  hurry  the 
operation  at  this  point. 

(5) .  After  all  the  ammonia  has  fumed  out,  the  current  is  further 
increased  carefully  to  melt  the  PbS  onto  the  filament.  Let  the  filament 
stay  at  a  red-hot  temperature  for  about  5  seconds  being  careful  not  to 
burn  out  the  filament. 

3»5«  Blanks 

Three  blank  runs  were  made  using  Pb^^-  as  a  spike.  They  were 

made  with  all  components  set  up  in  operating  condition  without  the 

sample  in  the  graphite  crucible.  The  furnace  was  heated  to  1400°C 

for  two  hours  and  the  deposit  on  the  cold  finger  was  dissolved  in 

on  L 

5  c.c.  of  3N  HNO^.  A  known  amount  of  Pb  spike  of  known  composition 
was  added  to  the  solution  and  the  lead  was  extracted  using  the 
dithizone  technique  with  the  same  amount  of  chemicals  as  normally 
used  in  processing  a  regular  sample.  Results  of  the  blank  tests  are 
listed  in  Table  3» 

It  can  be  seen  from  Table  3  that  the  blanks  are  approximately 
0.2  ^tg,  most  of  the  Pb  contamination  probably  coming  from  the  reagents 
used.  The  sample  size  for  most  of  the  runs  varied  from  over  50yng  to 


. 

' 


!  Li  ;  if  (  5  \:f  lo  teom  M  s. 


Table  3 •  Blank  tests  made  on  combined  volatilization 
and  dithizone  purification  of  lead. 
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to  about  20  fJ-g  estimated  from  the  amount  of  dithizone  used  for  the 
extraction.  No  blank  corrections  were  made  for  ratios  obtained  in 
this  thesis  as  it  was  felt  no  significant  change  in  the  interpreta¬ 
tion  of  the  results  would  arise  due  to  the  fact  that  the  isotope 
ratios  observed  lie  along  lines  which  pass  very  near  to  the 
composition  of  modern  lead  and  it  is  believed  that  the  blank  composi¬ 
tion  is  essentially  modern.  Thus  no  significant  change  in  the  slope 
of  the  line  through  the  various  isotope  ratios  will  be  produced. 

Separation  of  lead  by  volatilization  using  the  induction  furnace 
followed  by  dithizone  extraction  is  a  rapid  way  of  extracting  trace 
amounts  of  lead  from  rocks  or  other  material.  Lead  can  be  extracted 
and  ready  for  mass  spectrometric  analysis  within  a  few  hours.  If 
sampling  is  not  important,  the  rock  sample  can  be  heated  as  chips 
instead  of  powder  which  allows  the  evacuation  and  the  heating  to  be 
done  much  more  rapidly.  The  rock  chips  with  thickness  of  about  (1/4)" 
usually  were  completely  melted  within  15  minutes  at  a  temperature  of 
1300°C. 

All  the  leads  studied  in  this  thesis  except  those  from  galenas 
are  obtained  by  volatilization  followed  by  dithizone  extraction. 
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CHAPTER  IV 


MASS  DISCRIMINATION  IN  A  MASS  SPECTROMETER 


4.1.  Introduction 

The  basic  function  of  the  mass  spectrometer  is  to  produce  ions 
from  a  sample,  then  separate  them  according  to  their  mass  to  charge 
ratio.  All  practical  mass  spectrometers  used  to  determine  lead  isotope 
abundances  are  fundamentally  similar  (Russell  and  Farquhar,  i960). 

They  differ  only  in  the  methods  employed  in  producing  ions,  in  the 
method  of  measuring  and  recording  the  electron  current  which  neutra¬ 
lizes  the  ion  beams,  and  in  the  electronic  circuits  which  are  used  to 
produce  and  control  the  various  potentials  and  magnetic  field  by 
which  analysis  of  the  samples  is  accomplished. 

Most  investigators  concede  that  in  any  mass  spectrometer  there 
are  a  number  of  sources  of  mass  fractionation  which  might  introduce 
errors  into  the  final  results  (Russell  and  Farquhar,  i960).  In  a  gas- 
source  mass  spectrometer,  isotopic  fractionation  is  due  to  molecular 
flow  from  the  source  ionizing  region  to  the  main  volume  of  the 
instrument.  This  effect  should  be  constant  and  according  to  Whittles 
(1964)  for  samples  in  the  mass  range  of  lead  tetramethyl  it  amounts  to 
about  0.2  %  per  unit-mass  (quoted  from  Stacey  et  al,  1969).  For  a 
solid-source  mass  spectrometer,  various  sources  of  isotopic  fractiona¬ 
tion  have  been  suggested.  Cooper  and  Richards  (1966)  found  that  the 
beam-building  processes  can  often  have  a  significant  effect  on  the 
measured  ratios  (i.e.  fractionation  occurs  during  evaporation). 

Inghram  and  Hayden  (1954)  investigated  the  response  of  electron 
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multipliers  to  ion  mass  and  ion  velocities  and  according  to  them  the 
response  is  proportional  to  the  square-root  of  mass.  Doe  et  al  (1965) 
concluded  that  mass  dependent  variations  in  the  isotope  ratios  due  to 
different  materials  of  the  source  filament  were  observed  and  addition¬ 
al  mass  dependent  variations  exist  which  may  be  due  to  nonuniform 
sample  purity,  particularly  the  presence  of  alkalies.  On  the  other 
hand  Cooper  and  Richards  (1966)  failed  to  find  any  significant  effect 
from  different  chemical  procedures  on  the  observed  ratios.  In  their 
paper  on  triple- filament  lead  isotope  measurements,  Stacey,  Delevaux 
and  Ulrych  (1969)  found  that  the  exact  values  of  the  ratios  recorded 
in  any  one  analysis  seem  to  be  dependent  on  the  type  of  filament 
material  used  and  indirectly  on  the  exact  relative  alignment  of  the 
filaments  (triple- filament)  themselves. 

Readers  are  reminded  that  the  mass  discrimination  error  discussed 
in  this  thesis  refers  to  the  mass  dependent  errors  between  runs  on 
a  sample.  No  noticeable  mass  discrimination  between  sets  of  data 
taken  within  a  run  was  ever  observed  by  the  author  in  this  investi¬ 
gation.  Doe  et  al  (1965)  reported  that  no  measureable  changes  in  ratios 
were  noted  in  several  runs  for  which  two  or  three  sets  of  data  were 
taken  at  hourly  intervals.  Stacey  et  al  (1969)  stated  that  provided 
analyses  are  made  under  identical  conditions  using  the  same  filament 
material,  the  long-term  reproducibility  (i.e.,  over  a  period  of  months) 
is  not  appreciably  worse  than  the  day  to  day  reproducibility.  Cooper 
and  Richards (1966)  also  observed  that  no  appreciable  mass  discrimina¬ 
tion  occurred  during  the  data  taking  period  .  Catanzaro  et  al  (1968) 
noted  that  a  email  amount  of  isotope  fractionation  generally  occurs 


during  an  analysis.  According  to  them  the  Pb2^/Pb2<^  ratio  generally 
increased  by  about  0.03 %  during  the  data  taking  period  of  24  minutes. 
This  author  believes  the  fractionation  observed  by  Catanzaro  et  al 
might  have  resulted  from  a  slight  change  in  filament  conditions  which 
in  turn  altered  the  path  of  the  ion  beam  resulting  in  a  change  in  the 
observed  ratios. 

4.2.  Measuring  error  and  fractionation  due  to  the  positioning  of 

the  filament  in  the  mass  spectrometer 

In  order  to  obtain  estimates  of  the  measuring  error,  a  sample 
was  run  continuously  for  six  hours  obtaining  200  scans  of  the  lead 
mass  spectrum.  The  set  of  200  scans  was  subdivided  into  20  groups 
of  10  scans  each  and  ratios  were  calculated  using  5th  order  of 
polynomial  fitting  as' described  in  Chapter  II.  Results  are  shown  in 
Table  4* 

It  can  be  seen  that  the  scatter  is  quite  small,  about  0.06$  at 
one  standard  deviation  and  virtually  constant  for  each  of  the  three 
ratios  Pb  /Pb  ,  Pb  /Pb  and  Pb  /Pb  .  There  is  no 
observable  difference  between  the  ratios  from  the  first  group  and 
those  from  the  last;  that  is,  no  measureable  fractionation  has  taken 
place  for  this  particular  sample  during  the  six  hour  period  within  the 
above  mentioned  limits. 

The  above  sample  was  removed  from  the  mass  spectrometer  and  then 
replaced  several  days  later.  At  this  time  it  was  again  found  that  the 

ratios  obtained  were  constant  for  many  hours,  but  the  ratios  were 


higher  than  the  previous  long-term  run.  Results  of  this  run  are  listed 
in  Table  5.  This  suggests  that  the  positioning  of  the  filament  is 
very  critical.  To  further  test  this  possibility  the  above  mentioned 
sample  was  left  under  vacuum  in  the  mass  spectrometer  overnight  and 
then  the  filament  reheated  and  the  ratios  remeasured  the  next  day. 

It  was  found  that  the  ratios  were  the  same  as  those  shown  in  Table  5. 
Again  no  time  dependent  fractionation  was  observed.  Results  of  this 
measurement  are  shown  in  Table  6.  Thus  the  only  remaining  variable 
factor  seems  to  be  the  filament  position. 

Within  a  run,  the  position  of  the  filament  is  fixed  provided 
the  filament  does  not  alter  its  shape  during  the  run.  In  this  case 
the  measuring  error  is  independent  of  the  isotopic  mass  as  demon¬ 
strated  by  the  data  shown  in  Tables  4,  5  and  6.  Between  runs  the 
filament  positions  might  not  be  the  same  and  this  seems  to  introduce 

mass  dependent  errors  in  the  measured  ratios  as  shown  by  the  data  of 
Tables  1  and  2. 

It  is  known  that  gas-source  mass  spectrometry  has  better  repro¬ 
ducibility.  This  is  understandable  because  there  is  no  sample  filament 
involved  and  the  relative  position  of  the  source  to  collector  is 
always  the  same  between  different  runs,  which  is  in  fact  similar  to 
the  situation  for  a  set  of  data  obtained  within  a  run  of  a  single 
filament  solid— source  mass  spectrometer. 

4.3 >  Focusing  effects  on  the  isotonic  ratios 


To  see  if  focusing  has  any  effect  on  the  isotopic  ratios,  the 


Table  4*  Ratios  of  a  sample  continously  run  for  6  hours 

(Broken  Hill  #1  galena) 


,2°6  /pb20U 

Pb207/Pb20^ 

Pb208/Pb204 

16.066 

15.459 

35.844 

16.069 

15.480 

35.914 

16.070 

15.473 

35.918 

16.060 

15.464 

35.895 

16.061 

15.473 

35.926 

16.050 

15.462 

35.894 

16.045 

15.450 

35.876 

16.074 

15.474 

35.914 

16.074 

15.467 

35.883 

16.056 

15.467 

35.896 

16.059 

15.460 

35.878 

16.048 

15.454 

35.888 

16.060 

15.462 

35.836 

16.056 

15.460 

35.881 

16.071 

15.464 

35.869 

16.062 

15.474 

35.902 

16.074 

15.475 

35.902 

16.053 

15.457 

35.906 

16.073 

15.469 

35.396 

I6.O63 

15.467 

35.900 

Average:  16.062  +  .009  15.466  ±  .008  35.890  ±  .022 

(.056$)  (.045$)  (.061$) 
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Table  5.  Ratios  of  Broken  Hill  galena  (same  sample  as  that 
used  in  Table  4). 


Average 


Pb206/Pb2°4  Pb2Q7/Pb2°4  Pb2Q8/Pb2Q4 


16.139 

16.085 

16.115 

16.100 

16.085 

16.117 
16.134 
16.103 

16. 118 
16.087 
16.121 

16.109  +  .018 

(.11® 


15.555 

15.507 

15.528 

15.514 

15.512 

15.541 

15.547 

15.520 

15.534 

15.514 

15.541 

15.528  ±  .016 

(.11® 


36.122 

36.016 

36.086 

36.023 

36.042 

36.102 

36.080 

36.049 

36.068 

36.030 

36.084 

36.064  ±  .033 

(.09® 


Table  6.  Ratios  of  Broken  Hill  #1  galena  (same  sample  loading 
and  running  condition  as  that  of  Table  5) 


Pb206/Pb204 

Pb207/Pb2(* 

Pb208/Pb204 

16.126 

15.540 

36.112 

16.119 

15.529 

36.071 

16.090 

15.504 

36.013 

16.090 

15.503 

35.991 

16.119 

15.528 

36.073 

16.104 

15.530 

36.057 

16.087 

15.510 

36.044 

16.109 

15.524 

36.052 

16.112 

15.525 

36.058 

16.102 

15.520 

36.028 

16.089 

15.513 

36.015 

16.101 

15.526 

36.070 

16.116 

15.526 

36.038 

16.104 

15.522 

36.024 

16.111 

15.517 

36.020 

16.122 

15.537 

36.077 

16.108 

15.526 

36.077 

16.126 

15.542 

36.087 

16.101 

15.522 

36.059 

16.113 

15.529 

36.047 

16.098 

15.519 

36.053 

16.105 

15.506 

35.991 

Average:  16.107  ±  .011  15.523  ±  .010  36. 048  ±  .030 

(.062$)  (.062#)  (.0959 
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author  ran  one  Pb  spiked  sample  continuously  for  six  hours 
obtaining  200  scans  which  consist  of  20  sets  of  10  scans  each.  Among 
the  20  sets  of  data,  focusing  adjustments  were  made  on  Pb208,  Pb20^, 

206  202l 

Pb  and  Pb  for  each  5  sets  of  data  respectively.  The  sample  was 
spiked  with  Pb20^  to  increase  the  Pb20^  signal  so  that  the  Pb20^ 


measurement  errors  are  comparable  to  the  errors  on  the  other  isotopes. 
Results  of  this  test  are  shown  in  Table  7  and  shows  that  the  ratio 
pb2°8/pb204  observed  when  the  beam  was  focused  on  pb20S  ds  bhe  highest 
when  compared  to  the  same  ratios  when  the  beam  was  focused  on  other 
isotopes.  A  similar  result  was  also  noted  for  Pb20^  focusing.  The 
results  also  show  that  when  the  beam  was  focused  on  Pb2^^,  the  ratios 
Pb208/Pb204,  Pb207/Pb204  and  Pb206/Pb204  are  lower  when  compared  to 
the  same  ratios  observed  when  the  beam  were  focused  on  other  isotopes. 

Mass  spectrometric  analyses  using  a  triple-filament  (Stacey  et 
al,  1969)  and  double  spiking  technique  (Compston  and  Oversby,  1969) 
are  aimed  at  improving  the  reproducibility  of  mass  spectrometer 
measurements.  The  triple  filament  with  a  large  load  of  sample  and 
lower  filament  current  may  have  the  effect  of  making  the  path  of  the 
beam  more  constant  than  single  filament  measurements  and  consquently 
yield  better  reproducibility. 

In  the  double  spiking  technique,  the  mass  discrimination  effects 
are  corrected  by  making  an  additional  run  of  the  sample  mixed  with 
a  spike  of  known  composition  but  in  the  author’s  opinion  both  runs 
should  be  made  by  focusing  on  the  same  isotope  since  this  type  of 
discrimination  may  be  non-linear  being  affected  by  the  beam  focusing 


as  shown  above. 


fSO  i 
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4«4«  Conclusions  on  the  mass  discrimination  experiment 

Based  on  limited  results  and  discussions  in  this  Chapter,  the 
following  conclusions  are  made; 

(1) .  The  mass  discrimination  error  is  the  basic  limitation  on  the 

overall  reproducibility  of  a  mass  spectrometer, 

(2) .  The  mass  discrimination  error  is  related  to  filament  position 

in  the  mass  spectrometer  and  better  reproducibility  could  be 
obtained  if  a  means  of  making  the  sample  position  constant 
could  be  found, 

(3) .  No  noticeable  fractionation  during  the  course  of  a  six  hour 

run  was  observed  above  a  limit  of  about  0.06$  (S.D.).  Thus 
during  the  course  of  a  normal  run  lasting  one  or  two  hours, 
any  fractionation  would  be  negligible. 
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CHAPTER  V 


WHOLE  ROCK  DATA  AND  U/Pb  AGES 

5*1.  Whole  rock  Pbt  U  and  Th  determinations 

All  the  whole  rock  leads  except  that  of  YKWG  #8  were  extracted  by- 
heating  about  10  grams  of  sample  in  the  induction  furnace.  The  tempera¬ 
ture  was  raised  slowly  to  12+00° C  and  maintained  at  that  temperature 
for  approximately  30  minutes.  Appreciable  amounts  of  lead  are  observed 
as  soon  as  the  temperature  exceeds  about  900°C  and  hence  it  is  felt 
that  it  is  not  necessary  to  maintain  the  temperature  for  more  than  30 
minutes  at  12+00°C  in  order  to  obtain  complete  recovery  of  the  lead. 

Lead  collected  on  the  cold  finger  was  dissolved  in  10  c.c.  of 
3N  HNO^  and  split  into  two  portions  of  5c. c.  each,  one  for  lead-isotope 
composition  measurement  and  one  for  isotope-dilution  determination  of 
lead  content.  Purification  of  lead  was  done  by  the  dithizcne  extraction 
technique  as  described  in  Appendix  III. 

For  sample  YKWG  #8,  two  separate  heatings  one  with  plain  rock 
powder  for  lead  composition  and  one  with  Pb  4  spiked  powder  for  lead 
concentration  were  made. 

Replicate  volatilization  of  the  sample  YKWG  #6  was  made  to  give 
an  indication  of  the  furnace  reproducibility. 

U  and  Th  determinations  were  made  following  the  chemical  procedures 
described  in  Chapter  II.  Neutron  activation  analyses  were  also  made 
for  U  and  Th  determinations  (by  Dr.  Gumming). 

Tables  8  and  9  show  the  analytical  data  for  Pb,  U  and  Th  deter¬ 
minations.  Both  the  observed  and  spiked  lead  isotope  ratios  are  given. 
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Table  8.  Whole  rock  isotopic  ratios 

Pb  ratios,  unspiked  Pb  ratios,  204  spiked 

Sarngle:  206/204  207/204  208/204  207/206  208/206  206/204  207/204  207/206  208/206 
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Ratios  were  reduced  to  conform  to  NBS  Equal  Atom  standard  given  by 
Stacey  et  al  (1969). 

Table  10  shows  a  comparison  of  U  and  Th  concentration  determina¬ 
tions  by  the  isotope  dilution  measurements  and  neutron  activation 
analyses  on  the  same  rock  samples.  The  agreement  is  not  good  and  the 
reasons  for  the  disagreement  are  not  known.  The  isotope  dilution 
measurements  have  been  used  in  this  work  for  all  the  samples  except 
YKWG  #6,  when  the  isotope  dilution  measurement  was  not  available. 


5.2. 


Pb 


208 


Pb 


206 


U 


238 


Pb 


206 


Uranium-lead  radiometric  method  (after  Kanasewich,  1968) 

2o l  206  20V 

There  are  four  stable  isotopes  of  lead,  Pb  1  Pb  ,  Pb  and 
on  L 

.  Pb  4  has  not  been  generated  by  any  radioactive  decay  while 
207  208 

,  Pb  '  and  Pb  are  being  added  as  they  are  the  end  products  of 
U235  and  Th232  decays.  Therefore  the  relative  abundances  of 
,  Pb2^  and  Pb2<33  with  respect  to  Pb2<3^  will  increase  as  time 


progresses. 

Let  the  present  be  t=0  and  count  time  positively  into  the  past, 

and  let  N  be  the  number  of  radioactive  atoms  at  the  present  time, 
o 

Then  the  number  of  the  atom  in  the  past  at  time  t  is  N=NQe  ,  where 
X  is  the  decay  constant. 

With  respect  to  Pb2<3\  let 

a  ,  b  and  c  be  the  ratios  of  Pb2^/Pb2(3^,  Pb2<3^  and 
o’  0  o  ' 

pb208/pb204  reSpectively  at  t=tQ  (taken 
to  be  the  age  of  the  earth) 

?  Vn  238  235 

X ,  X  and  ^  be  the  decay  constants  of  U  ,  U  and 


Table  10.  Concentrations  of  U  and  Th 


U(ppm) 


Th(ppm) 


Sample : 

* 

I.  D. 

@ 

D.N. 

* 

I.D. 

@ 

D.N. 

YKWG  #1 

— 

7.62 

.14 

— 

20.0 

1.6 

YKWG  #2 

2.10 

2.08 

.04 

9.41 

9.36 

.70 

YKWG  #3 

3.19 

4»34 

.03 

15.98 

9.89 

.11 

YKWG  #4 

0.734 

0.62 

.01 

5.86 

5.43 

.02 

YKWG  #5 

1.20 

1.33 

.02 

5.80 

8.93 

.41 

YKWG  #6 

— 

13.51 

.15 

— 

40.9 

.3 

YKWG  #7 

1.49 

2.63 

.01 

3.76 

5.82 

.24 

YKWG  #8 

2.11 

2.33 

.01 

14.33 

12.45 

1.19 

*  Isotope  dilution  measurement  by  the  author. 

@  Delayed  neutron  measurement  by  Dr.  Cumming  (personal 
communication).  The  values  shown  are  the  means  of 
two  measurements.  The  errors  are  standard  deviation 
of  the  mean. 
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Th232  respectively, 

p.{t)  and  w(t)  be  the  ratios  of  the  number  of  atoms  of 
U238  and  Th232  with  respect  to  Pb204  at 
the  present  time,  and 

x,  y,  z,  be  the  ratios  of  Pb206/Pb2°4,  Pb207/Pb2°4  and 
Pb208/Pb204  at  time  t. 


Then  the  evolution  of  lead  since  the  earth  was  formed  can  be  expressed 


as 


x  =  a 


-p 

1 — 1 
-p 

-p 

i _ i 

3 

(4) 

t 

(1/137.8)  /  °^[ji(t)ex,t]dt  - 

- (5) 

t 

J  A”[w(t)eA  ]dt 
"t 

- (6) 

The  parameters  u(t)  and  w(t)  may  vary  from  point  to  point  in  the 
earth  but  metamorphic  or  orogenic  events  are  episodic  locally  and 
brief  in  comparison  to  the  half-lives  of  uranium  and  thorium,  so 
they  can  be  treated  as  constant  for  discrete  intervals  of  time  and 
approximate  solutions  to  the  above  equations  (4),  (5)t  and  (6)  can 
be  written  as 


/  At 

x  =  a  +  li- (e  o 
n  o  ^1 


>\t. 


e  l)  +  ^(e^l  -  eAU2)  +  .... 
+  ....+  nn(e  n-1  -  e  n) 


xt. 


(7) 


t  'tc  J  SC  ( 


■ 


n  s  —  I-n  si)  ^4  +■••«*  + 


yn  =  bo  +  (1/137.8)[  n^e 


•  • 


z 

n 


c 

o 


+  wx(e 


x"t 

o 


x’t  x't.,*  /  ^'t1  >’t0\ 

o  -  e  1)  +  p2(e  1  -  e  2)  + 

+  ...  +  pn(ex  ^n-1  -  e^^n)  ] - (8) 

eV,tl)  +  „2(e^"tl  -  6^2)  +  ... 

+  ...  +  wn(eX  ^n-1  -  ex*n)  - (9) 


where  n  is  the  number  of  stages  where  the  uranium  and  thorium  in  the 

system  has  been  altered  since  the  earth  was  formed  at  t  . 

o 

(A).  For  a  single  stage  lead 


X1  =  ao  +  jVeXt°  "  e>kt - (10) 

71  =  bQ  +  (l/l37.8)[ji1(ex,to  -  eX,tl) - (ll) 

Z1  =  co  +  wi^eX  to  ~  - (12) 

The  above  equations  may  be  used  to  solve  for  t^,  the  time  when  the 
system  became  closed. 

*1  -  b0  n  (e^o-e^l)  _ (13) 

X1  “  ao  137.8(ext0  -  extl) 

This  equation  defines  a  straight  line  with  slope  R  and  passes  through 
the  points  (aQ,  b  )  and  (x^,  y^).  All  leads  which  have  evolved  from 
the  same  initial  ratios  (Primordial  lead  )  through  different  environment 
which  became  closed  systems  at  the  same  time  t^  will  fit  on  the  same 
straight  line  called  an  isochron. 


XXxw 

•  IKnrfpC?  :  !'.*  i  S-  '  J' ; 
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(B)  .  For  two  stage  lead 

At  time  t^  the  leads  were  differentiated  from  their  source  and 
contemporaneously,  uranium  and  thorium  were  also  incorporated  into 
the  same  enviroment.  At  time  t^  a  geological  event  occurred  in  this 
area.  Some  of  the  initial  lead  formed  from  t  to  t.  and  the  radiogenic 
lead  formed  between  time  t^  and  t^  may  become  mixed.  Leads  which 
evolved  from  the  initial  ratio  (x^,  y^)  through  different  enviroments 
but  entered  closed  systems  at  the  same  time  t^  will  fit  on  a  straight 
line  defined  by 

y2  -  yl  (e^l-6^2)  _ (14) 

x2  X1  137. 8(e  tl  -e  t2) 

called  an  anomalous  lead  line. 

The  two  stage  model  assumes  that  the  incorporations  of  lead, 
uranium  and  thorium  into  the  crust  occurred  contemporaneously  during 
an  early  period  of  tectonic  activity  at  time  t^.  A  schematic  diagram 
illustrating  the  development  of  lead  and  a  lead-lead  isochron  diagram 
of  the  two  stage  model  are  shown  in  Figure  11. 

(C) .  Three  stage  lead  evolution 

For  three  stage  lead,  we  have 

x^  =  aQ  +  p^e^o  -  e^l)  +  p^e^l  -  e  X*2) 

+  p^(ext2  -  ext3)  - (15) 

y3  =  bQ  +  (1/137.8)[  p^e^o  -eVtl)  +  p2(eVtl  -  ex,t2) 

+  p3(ex  t2  -  e*  t2)]  - (l6) 


P21 

P22 

^23 

Pi 

• 

e 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

*2n 

Event :  Primordial  Pb 

and  U  introduced 
into  system 


Episodic  change 
in  U/Pb  ratio 


U,  Pb  analyses  of 
samples  at  the 
present  day 


Figure  11.  The  development  of  lead  in  a  two  stage  system  and 

the  lead-lead  diagram  (after  Gale  and  Mussett,  1973)* 
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z3  =  co  +  w^e^o  -  e^l)  +  w2(ex''Tjl  - 


\"t. 


x”t 


K"t 


2) 


+  w  (^2  -  e^3) 


•(17) 


If  the  uranium  and  thorium  concentration  varies  from  place  to  place, 
the  isotopic  ratios  will  scatter  about  an  anomalous  lead  line  with 
the  slope  determined  by  times  t^  and  t^. 

Only  under  two  special  cases  when  anc*  ^2i/^3j  ^or  eac^ 

of  the  subsystems  are  constants  (see  Figure  12)  will  a  rectilinear 
plot  of  Pb  /Pb  ^  versus  Pb^^/Pb^^  be  obtained  (Gale  and  Mussett, 
1973).  In  general  it  is  impossible  to  obtain  a  unique  solution  of  t  , 
^1*  ^2*  ^1  an<^  ^2  ^or  three  stage  lead  from  the  isotope  data  alone. 

For  the  two  special  cases,  Gale  and  Mussett  (1973)  have  shown 
that  for  the  case  fli=M  l/^2i=C0ns'tant ,  an  isochron  equation 


where 


y3  -  y" 


(ex,t2  -  ex,t3) 


137.8(eXt2  -  exl,3) 


■(18) 


xjb. 


xt 


At. 


x"  =  a0  +  ^(a"bo  -  e^l)  +  (p1/f1)(ext’l  -  eXt2) 


At. 


y”  =  bQ  +  (1/137.8)[  ^(eN,to  -  e A' X 1 )  +  (n/f ^(e^l  -  eVt2)] 


X't. 


x’t. 


Vt, 


is  obtained.  These  isochrons  define  a  family  of  parallel  lines  with 
meaningless  intercepts  with  the  growth  curve.  The  slope  of  the  isochron 
depends  on  t^  and  t^  but  not  on  t^.  When  X,  is  assigned,  t^  is  calcu¬ 
lated  from  the  slope.  As  f^  approches  1,  the  model  reduces  to  a  two 
stage  model.  Figure  12  shows  the  schematic  diagram  illustrating  the 
development  of  this  particular  model  lead  and  the  associated  isochron 
diagram. 


'■ 

. 


•-,i : v.  ,  •  ■  rt  •  -  1  - 


Figure  12.  The  development  of  lead  in  a  three  stage  system  and  the 
lead-lead  diagram  for  lead  that  has  evolved  according 
to  a  three  stage  f -^constant  history. 

(atfer  Gale  and  Mussett,  1973) 


. 
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For  the  case  f, 


^constant ,  the  isochron  equation  is 


given  by 


where 


=  ao  +  u1(eXto  -  extl) 
y’  =  bQ  +  (u1/l37.8)(e>''to  -  eVtl) 


In  this  case  the  possible  isochrons  have  different  slopes  (dependent 
on  the  value  of  f^,  but  all  pass  through  the  point  (x^ty^)  which  is 
the  point  of  intersection  of  the  single  stage  growth  curve  with  the 
single  stage  t^  isochron.  The  lead-lead  diagram  of  this  particular 
case  of  lead  evolution  is  shown  in  Figure  13 . 


(D) .  Uranium-lead  and  thorium-lead  ages 


Determination  of  uranium-lead  and  thorium-lead  ages  of  a  rock 
depends  on  a  determination  of  the  total  amount  of  radiogenic  lead 
produced  by  the  decay  of  the  parent  U  or  Th  subsequent  to  an  initial 
state  of  the  rock  assuming  it  to  be  a  closed  system  thereafter.  Three 
independent  ages  are  obtained  from  the  following  three  equations 


(20) 


-(21) 


'  . 


. 


-  . 


(  ) 


. 
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Figure  13.  Lead-lead  diagram  for  lead  that  has  evolved  according  to 

a  three  stage  f2=constant  history.  The  region  of  the 

diagram  accessible  to  the  tertiary  isochrons  is  limited 

in  both  the  f 2=Qo  and  f  2=0  directions . 

(after  Gale  and  Mussett,  1973) 
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'208 


-  OK208  208>. 

= J_  log  (  1  +  a _ -  (Fb  )t-t.1  \ 

*  6  (Th232)t=0 


■(22) 


206 

where  Pb  is  the  amount  of  lead  206  present  in  the  rock  at  the  time 
of  analysis  and  (Pb  )t=ti  the  amount  incorporated  into  the  rock 
when  the  system  became  closed.  For  a  two  stage  evolution  of  lead  the 
plot  of  Pb^  versus  Pb^^/U^  ^  is  a  straight  line  which  inter¬ 

sects  the  concordia  curve  at  two  points  (Wetherill,  1956;  Ulrych, 

1967)  as  shown  in  Figure  14  .  The  lower  intercept  yields  t^,  the  time 
when  the  system  became  closed  and  the  upper  intercept  yields  an 
independent  age  of  the  earth.  In  general  individual  Pb'^/U^^  and 
Pb^ytn^  ages  are  usually  discordant  due  to  the  loss  of  lead 
and/or  loss  or  gain  of  uranium.  Wetherill  (1956)  showed  that  if  Pb 
loss  occurs  in  a  single  episode  after  the  original  formation  of  the 

U-bearing  mineral,  then  a  concordia  plot  of  Pb'^/U*^  versus 
206  238 

Pb  /U  for  cogenetic  samples  of  varying  Pb  loss  will  yield  a 
straight  line  which  intersects  the  concordia  curve  at  the  time  of 
the  original  mineral  formation  and  the  time  of  the  subsequent  episodic 
lead  loss  (Baadsgaard,  1965). 

According  to  Tilton  (i960)  Pb  may  diffuse  slowly  from  a  mineral 
over  its  lifetime,  rather  than  being  lost  at  one  or  more  discrete 
intervals.  A  plot  of  this  model  of  Pb-loss  on  a  concordia  diagram 
gives  a  straight  line  over  most  of  the  concordia  plot,  but  the  line 
curves  to  the  origin  at  low  values  of  Pb^6/u^38  p^207^y235 

the  upper  intercept  being  the  time  when  the  system  became  closed. 


206  /tt238 
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o 

UA 


Figure  1 4.  Modified  concordia  plot  for  a  two-stage  evolution 
of  lead,  (after  Gale  and  Mussett,  1973) 


For  the  two  special  cases  of  three-stage  lead,  f ^constant  and 
f2=constant,  Gale  and  Mussett  (1973)  showed  that  the  modified  con- 
cordia  plots  are  families  of  straight  lines  as  shown  in  Figures  15a 
and  15b.  In  the  case  of  f  = const ant,  all  the  lines  pass  through  the 
same  lower  intercept  with  the  concordia  curve  at  t2  but  in  the  case 
of  f2=constant,  all  the  lines  are  parallel  with  meaningless  upper 
and  lower  intercepts,  and  the  slope  of  the  line  yields  t^  if  tQ,  the 
age  of  the  earth,  is  assumed. 

5,3.  Whole  rock  lead-lead  and  uranium-lead  ages 

Figure  16  shows  the  plot  of  Pb2^/Pb2^  versus  Pb  ^/Pb  ^  for 
the  whole  rock  lead  isotope  data  of  Table  8.  The  best  straight  line 
was  calculated  following  the  method  of  Cumming  et  al  (1972)  using 
equal  weighting  factors  for  all  ratios.  The  slope  of  this  line  as 
given  by  equation  (14)  is 

R  =  0.1808  +  .0094  (one  standard  deviation). 

The  two-stage  lead  model  age  at  which  the  system  had  an  initial  uni¬ 
form  lead  composition  is  determined  to  be  2700  ±  75  m.y.  (S.D.).  This 
age  is  slightly  higher  than  the  Rb/Sr  whole  rock  isochron  age  of 
2610  +  58  m.y.  (Green  et  al,  1968)  for  the  same  intrusive  rock.  The  p. 
value  for  the  growth  curve  which  intersects  the  above  line  at  t^=2700 
m.y.  and  t2=0  is  calculated  to  be  8.8  ±  .10  (95$  C.L.)  similar  to  the 
value  for  the  comformable  leads  obtained  by  Cooper  et  al  (1969). 

The  individual  Pb2^/U2"^  and  Pb  ^/U  ^>'>  ages  are  discordant 
(see  Table  9)  and  follow  a  general  discordant  trend  of  lead  loss  with 
the  Pb207/Pb206  isochron  age  being  higher  than  the  Pb20^/U2^  age 
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Figure  15a.  Modified  concordia  diagram  for  lead  that  has  evolved 
according  to  three- stage  f ]= const .  history. 

(after  Gale  and  Mussett,  1973) 


Figure  15b.  Modified  concordia  diagram  for  lead  that  has  evolved 
according  to  a  three-stage  f2=constant  history.  The 
accessible  region  of  the  diagram  is  limited  only  in 
the  f 2=0  direction,  (after  Gale  and  Mussett,  1973) 
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which  in  turn  is  higher  than  the  Pb  /U  age.  These  discordant  ages 
may  indicate  that  the  samples  have  lost  a  fraction  of  their  lead  con¬ 
tent  before  or  during  their  analysis  in  the  laboratory.  The  effect  o£ 
this  loss  is  to  lower  the  Pb  /U  age  proportionately  more  than 
Pb  °7/U235  aSe»  but  since  the  Pb2^7/Pb2^  age  is  based  entirely  on 
the  relative  proportions  of  two  isotopes,  their  ratio,  and  hence  the 
age  calculated  from  that  ratio  will  be  relatively  insensitive  to  loss 
of  lead  (Russell  and  Farquhar,  i960). 

The  modified  concordia  plot  of  the  form 

(Pb206/U238)o*=[(Pb206/Pb20bo  -  ao](Pb2CV38)0 

versus 

(Pb207/U235)*=[(Pb207/Pb204)o  -  bo](Pb204/U235)o  is  shown  in  Figure 
17. 

A  linear  trend  is  observed  although  the  data  show  considerable 
scatter,  and  no  simple  model  can  account  completely  for  the  data.  We 
may  however  fit  a  straight  line  to  the  data  points  (omitting  sample 
point  YKWG  #2)  and  this  line  intersects  the  normal  concordia  curve 
at  1540  m.y.  and  5000  m.y.  These  intercept  ages  cannot  be  interpreted 
on  the  assumption  of  an  approximation  to  a  two  stage  model  and  the 
upper  intercept  age  of  5000  m.y.  makes  it  impossible  to  apply  the 
single  episodic  lead-loss  theory  proposed  by  Wetherill  (1956)  or  the 
continous  loss  of  lead  by  diffusion  described  by  Tilton  (i960). 

Because  of  the  lack  of  a  consistent  pattern  in  ages  obtained  from 
different  approaches, that  is  the  lead-lead  isochron  diagram,  the  U-Pb 
age  equations  and  the  concordia  diagram,  it  is  clear  that  a  simple 
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two-stage  model  is  not  adequate  to  describe  the  lead  evolution  in  this 
area  even  as  a  first  approximation  although  the  lead-lead  isochron  age 
of  2700  m.y.  ±  85  m.y.  agrees  within  the  statistical  limits  of  the 
Rb/Sr  age  of  2610  m.y.  +  58  m.y.  obtained  by  Green  et  al  (1968).  It 
was  then  decided  to  perform  some  studies  of  fractional  removal  of  lead 

from  each  of  the  eight  sample  rocks  by  volatilization.  Cumming  et  al 

\  1 

(1970)  have  shown  that  incremental  heating  of  a  rock  sample  yields 
lead  of  varying  isotopic  composition  depending  on  the  heating  tempe¬ 
rature  and  the  data  fit  on  a  single  straight  line  analogous  to  a 
mineral  isochron.  It  was  hoped  that  this  inhomogeneous  distribution  of 
lead  in  rocks  would  help  to  clarify  the  age  differences  obtained  by  the 
different  methods.  The  results  of  this  volatilization  studies  are 
given  in  Chapter  VI. 

5.4.  Possible  three  or  higher-stage  system 

An  examination  of  the  concordia  plot  of  Figure  17  shows  that 
the  best  straight  line  which  fits  the  six  data  points  is  parallel  to 
one  which  has  an  upper  intercept-age  of  4578  m.y.  (the  age  of  the 
earth)  and  a  lower  intercept-age  of  2450  ±  450  m.y.  (S.D.)  with  the 
normal  concordia  curve.  If  a  three-stage  f ^constant  lead  evolution 
history  were  assumed  as  first  approximation  to  the  data,  the  lower 
intercept-age  calculated  from  the  slope  of  the  best  fitting  straight 
line  and  the  age  of  the  earth  is  t^  the  time  when  the  system  was 
subjected  to  the  first  change  in  U/Pb  ratios.  There  is  no  simple  way 
of  finding  t^  from  the  modified  concordia  diagram. 

It  is  not  clear  whether  a  three-stage  f2=constant  lead  evolution 
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model  can  be  used  to  represent  the  observed  data,  since  the  straight 
line  of  the  modified  concordia  diagram,  which  has  a  lower  intercept- 
age  of  1540  m.y.  +  200  m.y.,  lies  considerably  below  that  of  f£=l 
(see  Figure  17)  and  according  to  Gale  and  Mussett  (1973)  the  accessible 
region  is  limited  by  the  line  which  represents  f^O  anc^  has  a  tower 
intercept-age  of  t^.  Therefore  if  the  three-stage  f^  const  ant  model  is 
used,  the  time  t9  when  the  system  was  subjected  to  a  second  event 
which  resulted  in  the  changes  of  U/Pb  ratios  of  the  system  must  be 
less  than  the  above  lower  intercept-age  of  1540  m.y,  +  200  m.y..  The 
only  data  from  the  Yellowknife  area  which  may  fall  within  these  limits 
is  the  poorly  defined  Rb/Sr  isochron  on  dyke  set  III  (Gates  and  Hurley, 

1973 ) . 

5.5«  Second  whole  rock  ratios 

As  mentioned  in  section  5.1  the  lead  samples  were  extracted  by 
heating  rock  powder  at  1400°C  for  30  minutes.  For  the  purpose  of 
experimentation,  the  author  reheated  three  samples  two  of  these  in 
the  form  of  chips  instead  of  powder  (samples  YKWG  #3  and  #5)  and  one 
with  powder  to  check  the  reproducibility.  All  the  samples  were  heated 
at  1400°C  for  3  hours.  The  use  of  rock  chips  allows  faster  vacuum 
evacuation  and  faster  increase  of  heating  temperature.  All  chips  (up 
to  thick)  were  completely  melted  within  15  minutes  when  the  tempera¬ 
ture  rose  above  1300°C»  Table  11  shows  the  ratios  of  the  three  reheated 
samples  along  with  the  corresponding  ratios  from  Table  8. 

It  can  be  seen  that  ratios  from  the  two  separate  heatings  made 
on  powdered  sample  from  YKWG  #8  are  the  same.  The  differences  in  ratios 
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Table  11.  Comparison  of  whole  rock  Pb  ratios  between  the  first 
heating  of  Table  8  and  the  second  heating,  1973 


Pb206/^204 

Pb' 

207 /pb204 

Pb208/Pb2°4 

Sample : 

1973 

Table  8 

1973 

Table  8 

1973 

Table  8 

YKWG  #3 

16.336* 

21.445 

15.284* 

16.153 

36.407* 

36.395 

.017® 

.069 

.016 

.053 

.037 

.120 

YKWG  #5 

18.108* 

17.763 

15.612* 

15.450 

39.H2* 

38.052 

.024 

to 

<r\ 

O 

. 

.021 

.034 

.052 

.080 

YKWG  #8 

18.657 

18.654 

15.726 

15.748 

40.843 

40.291 

.032 

.103 

.027 

.088 

.072 

.228 

*  Sample 

was  heated 

as  rock  chips. 

@  Internal  estimate  of  standard  deviation 


for  samples  YKWG  #3  and  #5  are  believed  to  be  mainly  due  to  sampling 
problems  arising  from  the  fact  that  the  first  heating  was  made  with 
rock  powder  while  the  second  heating  was  made  with  chips. 

The  lead-lead  isochron  diagram  of  the  three  data  points  is  shown 
in  Figure  18.  A  sharply  defined  straight  line  with  slope  of 

R  =  0.188  ±  .008  (S.D. ) 

is  obtained.  The  slope-age  is  calculated  to  be  2760  m.y.  ±  85  m.y. 
which  agrees  within  the  statistical  limits  with  the  age  2700  m.y.  + 

85  m.y.  obtained  previously. 
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Figure  18 


CHAPTER  VI 


ISOTOPIC  COMPOSITION  OF  FRACTIONALLY  VOLATILIZED 


LEAD  AND  THE  TEMPERATURE  ISOCHRON 


6.1,  Introduction 

The  existence  of  different  modes  (forms)  of  lead  occurrence  in 
natural  environments  was  indicated  by  Tilton  (1956)  who  obtained  a 
displacement  in  the  isotopic  composition  of  lead  by  acid  leaching  of 
certain  minerals.  Starik  et  al  (1957)  conducted  extensive  volatili¬ 
zation  experiments  with  pitchblende,  monazite,  uraninite  and  a  series 
of  granites  and  made  the  following  conclusions. 

"/  \ 

(1) .  Lead  is  present  in  minerals  as  well  as  in  rocks  in 

different  forms. 

(2) .  Different  forms  of  lead  sublimate  with  different  speeds 

when  the  parent  minerals  or  rocks  are  heated.  This  results 
in  a  displacement  of  the  isotopic  composition  of  the 
isolated  lead. 

(3) .  Sublimation  in  vacuum  of  lead  from  pitchblende  (from 

Joachimsthal)  is  toward  an  impoverishment  of  common  lead 
in  proportion  to  an  increase  in  temperature  and  duration 
of  heating.  When  sublimated  in  hydrogen  the  shift  in  the 
isotopic  composition  is  toward  an  enrichment  in  radiogenic 
isotopes,  and  this  change  proceeds  much  more  smoothly. 
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(U) •  The  displacement  in  isotopic  composition  of  uraninite  is 
in  the  direction  of  enrichment  in  extraneous  lead,  namely, 
masses  204  and  208.  In  monazite,  on  the  other  hand,  the 
shift  is  toward  an  impoverishment  of  extraneous  lead. 

(5) «  When  a  sample  of  ancient  Karelian  granite  was  sublimated 

in  a  stream  of  hydrogen  and  the  Caledonian  granite  of 
Terek  Ala-Tau  in  vacuum, first  fractions  were  enriched  in 
radiogenic  isotopes.  Apparently,  radiogenic  lead  occurs  in 
granites  to  a  large  extent  disseminated  in  the  metallic 
state. 

(6) .  Under  natural  conditions  impoverichment  and  enrichment  of 

any  isotope  of  lead  regardless  of  its  mode  of  occurrence 
may  take  place  at  relatively  moderate  temperature  (about 
600°C). 

(7) .  Investigation  of  modes  of  lead  occurrence  by  sublimation 

methods  may  add  in  the  solution  of  problems  relating  to 
the  age  and  genesis  of  some  rocks  and  minerals." 

As  pointed  out  by  Starik  et  al  (1957),  the  different  origin  of 
lead,  determined  by  a  number  of  factors,  in  turn  determines  the 
various  modes  of  lead  occurrence.  The  modes  of  occurrence  of  primary 
lead  were  determined  by  physico-chemical  conditions  of  the  primeval 
material  that  make  up  the  earth.  As  a  result  of  natural  radioactive 
decay  of  uranium  and  thorium,  radiogenic  lead  is  being  continously 
generated  and  its  mode  of  occurrence  should  be  quite  distinct. 

Because  of  processes  of  remelting,  granitization,  weathering,  etc., 
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original  and  radiogenic  leads  may  occur  together  in  many  different 
ratios.  Identification  of  the  lead  types  should  determine  the  mode 
of  occurrence  of  common  lead.  The  chemical  compounds  in  which  common 
lead  participates  are  varied.  Radiogenic  lead,  on  the  other  hand, 
accumulating  in  a  rock  or  mineral  which  has  already  formed,  evidently 
remains  either  in  the  native  state  or  in  the  form  of  oxides.  Its 
position  in  the  crystal  lattice  of  a  mineral  depends  upon  the  crys- 
tallochemical  peculiarities  of  the  given  mineral.  Its  mode  of  occur¬ 
rence  depends  on  the  position  in  the  mineral  of  parent  elements 
related  to  its  original  formation  and  those  radioactive  elements 

related  to  secondary  mineralization.  However,  a  certain  portion  of 
206 

Pb  in  a  uranium  mineral  would  be  in  a  position  different  from  the 
bulk  of  radiogenic  lead  even  if  there  had  been  no  secondary  minerali- 

OO  I 

zation,  because  a  part  of  U  does  not  enter  the  crystal  lattice  ( 
Starik  et  al,  1957)*  The  same  may  be  said  of  some  other  intermediate 
decay  products  of  uranium  and  actino-uranium.  This  explains  the 
"fragile  structures"  of  radiogenic  lead  occurrence  (Starik  et  al, 

1957)- 

Generally  speaking,  processes  of  leaching,  emanation,  diffusion, 
volatilization,  recrystallization  and  isotopic  exchange  may  bring 
about  changes  in  the  isotopic  composition  of  lead  in  nature  and  inter¬ 
mixing  of  different  forms  of  lead. 

The  Volatilization  technique  has  been  used  for  some  time  as  a  means 
of  extracting  lead  from  minerals  and  rocks  (e.g.,  Edwards  and  Hess, 
1956;  Marshall  and  Hess,  I960;  Masuda,  1964;  Hamilton,  1965;  Cooper 
and  Richards,  1966;  Welke  et  al,  1968;  Cooper  and  Green,  1969;  Sinha, 


1969;  Cumming  et  al,  1970;  Black  et  al,  1971;  Silver,  1970  and  1971; 
Huey  et  al,  1971;  Doe  and  Tatsumoto,  1972),  but  up  to  the  present  time 
there  have  been  few  investigations  that  dealt  with  the  problem  of 
modes  of  lead  occurrence.  Different  modes  of  lead  occurrence  were 
indirectly  indicated  by  Sinha  (1969)  who  measured  isotopic  composi¬ 
tions  of  lead  in  three  fractions  of  potassium  feldspars  (the  bulk 
material,  the  first  volatile  fraction,  and  the  residual  volatile 
fraction)  from  granitic  rocks  and  pegmatites  and  found  that  lead  was 
more  radiogenic  in  the  first  volatile  fraction  and  less  radiogenic 
in  the  residual  volatile  fraction  compared  to  the  bulk  sample, 
indicating  that  radiogenic  lead  is  not  uniformly  distributed  in 
feldspar  crystals.  Sinha  interpretated  the  relative  ease  with  which 
the  radiogenic  lead  can  be  volatilized  as  a  suggestion  that  radio¬ 
genic  lead  is  concentrated  along  grain  boundaries  and/or  in  distorted 
lattice  sites  in  feldspar  crystals. 

Cooper  and  Green  (1969)  made  the  following  evaluation  of  the 
volatilization  separation  of  lead; 

If  the  efficiency  of  lead  volatilization  is  less  than  100^, 
it  remains  possible  that  the  data  obtained  are  not  represen¬ 
tative  of  the  whole  rock  lead  because  the  extraction  process 
has  preferentially  favoured  the  lead  in  some  minerals  relative 
to  others.  Such  partial  lead  extractions  should  yield  a  linear 
distribution  of  points  on  the  Pb^^/Pb^^  versus  Pb^^/Pb^^ 
diagram,  with  slope  dependent  on  the  age  of  crystallization 
of  the  rock  sample.  The  line  developed  would  extrapolate 


back  to  the  initial  lead  isotopic  ratios  of  the  parent  material." 

Cumming  et  al  (1970)  independently  observed  and  exploited  the 

inhomogeneous  distribution  of  lead  isotopic  composition  and  showed 

that  the  range  of  isotopic  data  on  a  plot  of  Pb^^/Pb^4'  versus 
207  204 

Fb  yPb  4  fit  on  a  simple  straight  line  from  which  some  information 
regarding  the  history  of  the  rock  may  be  obtained. 

Volatilised  lead  from  some  of  the  Apollo  12  and  14  soils  reported 
by  Doe  and  Tatsumoto  (1972)  showed  that  the  ratios  become  more  radio¬ 
genic  as  temperature  increases  in  the  range  600°C  to  1000°C  which  is 
contrary  to  the  results  observed  by  Cumming  et  al  (1970)  and  Black  et 
al  (1971)  for  terrestrial  rocks,  but  the  linear  relationship  between 
the  ratios  Pb^^/Pb^^  and  Pb^^/Pb^4-  is  evident,  except  for  one  data 
point  extracted  at  temperature  1350°C. 

6.2.  Temperature  isochron 

The  isotopic  compositions  of  leads  removed  from  the  eight  sample 
rocks  at  various  different  temperatures  are  shown  in  Table  12.  The 
leads  of  each  rock  sample  were  obtained  by  volatilizing  approximately 
20  grams  of  rock  powder  mixed  with  a  little  pure  graphite  powder  at 
each  of  the  various  indicated  temperatures  for  3  hours  in  the  high 
temperature  brick  furnace  (see  Chapter  III). 

It  can  be  seen  that  the  volatilized  leads  have  different  isoto¬ 
pic  composition  depending  on  the  heating  temperature  and  in  general 
lead  extracted  at  lower  temperature  is  more  radiogenic  than  that 
extracted  at  higher  temperatures.  There  is  also  evidence  that  when 
the  heating  temperature  rises  above  a  certain  critical  value,  which 
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Table  12.  Pb  ratios  of  fractionally  volatilized  lead 


Sample : 

Volatilization 
temperature  (°C) 

Pb206/Pb20^ 

Pb207/Pb20i 

pb208/pb204 

YKWG  #1 

750  -  850 

27.503 

.044* 

16.536 

.026 

40.256 

.064 

850  -  950 

24.222 

.033 

16.202 

.022 

38.702 

.053 

950  -  1050 

16.412 

.021 

15.232 

.020 

34.635 

.045 

1050  -  1150 

14.832 

.013 

15.135 

.013 

34.137 

.030 

Whole  rock 
(  1400° C) 

18.324 

.027 

15.627 

.024 

36.158 

.055 

YKWG  #2 

850  -  950 

20.593 

.094 

16.149 

.074 

41.998 

.198 

20.441 

.014 

15.980 

.011 

41.533 

.029 

950  -  1050 

18.015 

.015 

15.535 

.013 

38.023 

.031 

1050  -  1150 

16.537 

.017 

15.282 

.016 

36.630 

.038 

Whole  rock 
(  1400° C ) 

18.473 

.020 

15.817 

.017 

39.331 

.044 

*  Internal  estimate  of  standard  deviation. 


Table  12  coni’ d 


Sample : 


YKWG  #3 


YKWG  #4 


Volatilization 
temperature ( °C ) 

Pb206/Pb204 

Pb207/pb2C* 

Pb208/Pb20i4- 

750  -  850 

21.165 

.020* 

16.037 

.015 

37.203 

.036 

85 0  -  950 

17.636 

.016 

15.428 

.014 

35.277 

.032 

950  -  1050 

15.509 

.015 

15.060 

.014 

34.118 

.033 

1050  -  1150 

14.773 

.010 

14.978 

.on 

33.889 

.024 

Whole  rock 
( 1400° C) 

21. 446 
.070 

16.153 

.053 

36.395 

.119 

750  -  850 

20.110 

.868 

16.950 

.732 

51.789 

2.254 

850  -  950 

17.377 

.026 

15.390 

.023 

44.984 

.069 

950  -  1050 

16.498 

.029 

15.246 

.027 

39.780 

.072 

1050  -  1150 

16.868 

.027 

15.325 

.025 

36.127 

.061 

Whole  rock 
( 1400° C) 

16.504 

.102 

15.293 

.100 

40.265 

.265 

*  Internal  estimate  of  standard  deviation 
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Table 

12  cont'd 

Sample : 

Volatilization 
t emper atur e ( °C ) 

Pb^/Pb20^ 

Pb2°7/Pb20Zf 

Pb208/Pb: 

YKWG  #5 

750  -  850 

19.599 

.045* 

15.808 

.037 

41.716 

.099 

19-521 

.047 

15.758 

.038 

41.640 

.104 

850  -  950 

18. 514 
.025 

15.568 

.021 

39.927 

.055 

950  -  1050 

16.503 

.029 

15.255 

.027 

36.376 

.066 

1050  -  1150 

16.039 

.012 

15.187 

.011 

34.769 

.025 

Whole  rock 
( 1400° C) 

17.763 

.040 

15.450 

.035 

38.052 

.069 

YKWG  #6 

750  -  850 

60.878 

.051 

23. 088 
.020 

57.614 

.049 

950  -  1050 

28.049 

.047 

17.370 

.029 

40.736 

.069 

1050  -  1150 

20.891 

.029 

16.062 

.023 

36.724 

.052 

over  1150 

16.983 

.024 

15.489 

.022 

34.986 

.052 

Whole  rock 
(1400°C) 

40.802 

.140 

19.550 

.066 

45.876 

.157 

*  Internal  estimate  of  standard  deviation 


Table  12  cont’d 


Sample : 

Volatilization 

temperature(°C) 

Pb206/Pb204 

pb207/pb204 

Pb208/Pb2°4 

YKWG  #6 
(Reheat) 

750  -  850 

67.177 

.218* 

24.245 

.079 

58.422 

.193 

850  -  950 

49.782 

.099 

21.184 

.043 

50.116 

.102 

49.843 

.042 

21.241 

.018 

50.239 

.043 

950  -  1000 

42.155 

.023 

19.853 

.on 

46.528 

.026 

42.192 

.065 

19.867 

.031 

46.546 

.072 

1000  -  1150 

22.756 

.039 

16.459 

.028 

37.546 

.065 

22.798 

.007 

16.485 

.005 

37.607 

.012 

YKWG  #7 

750  -  850 

27.294 

.036 

16.941 

.023 

39.325 

.053 

850  -  950 

21.725 

.012 

16.100 

.009 

37.044 

.021 

950  -  1050 

18.070 

.017 

15.538 

.015 

35.491 

.034 

1050  -  1150 

16.877 

.011 

15.314 

.010 

34.966 

.023 

Over  1150 

15.520 

.024 

15.107 

.023 

34.162 

.053 

Whole  rock 
(1400°C) 

19.919 

.048 

15.968 

.039 

36.496 

.090 

Table  12  cont'd 
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Sample : 

Volatilization 

temperature(°C) 

Pb206/Pb204 

Pb207/Pb20^ 

Pb208/Pb204 

YKWG  #8 

750  -  850 

21.363 

.026* 

16.127 

.020 

51.006 

.062 

21.396 

.031 

16.136 

.024 

51.036 

.076 

850  -  950 

17.388 

.043 

15.338 

.038 

38.871 

.100 

950  -  1050 

16.632 

.046 

15.297 

.043 

— 

16.667 

.034 

15.407 

.032 

— 

16.743 

.037 

15.468 

.035 

— 

1050  -  1150 

23.751 

.119 

16.615 

.084 

— 

Whole  rock 
( 1400° C) 

18. 654 
.103 

15.748 

.088 

40.291 

.228 

*  Internal  estimate 

of  standard 

deviation. 

Oi.  ' 
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probably  depends  on  the  mineralogy  of  the  individual  sample  rock,  the 
volatilised  lead  becomes  more  radiogenic  again  at  the  higher  tempera¬ 
tures.  The  different  ratios  observed  for  different  extracting  tempera¬ 
tures  indicate  possible  differences  in  the  modes  of  lead  occurrence. 

The  fact  that  the  lowest  temperatures  yield  the  most  radiogenic 
lead  suggests  that  lead  generated  by  the  interstitial  material  which 
often  contains  appreciable  amounts  of  uranium  and  thorium  (Austin 
and  Slawson,  196 l)  breaks  down  first  and  subsequent  heatings  yield 
lead  from  the  more  refractory  minerals.  It  is  known  that  the  uranium¬ 
bearing  mineral  zircon  is  very  refractory  and  contains  highly  radio¬ 
genic  leads.  For  rocks  which  contain  abundant  zircon,  we  would  expect 
to  see  the  isotopic  composition  of  volatilised  lead  become  more 
radiogenic  as  the  temperature  exceeds  that  above  which  the  zircon 
starts  to  break  down.  This  phenomenon  is  evidented  in  samples  YKWG  #4 
and  #8. 

It  is  clear  that  lead  is  very  inhomogeneous  in  the  rock  as  a 
whole  and  that,  unless  extraction  is  quantitative,  the  volatilised 
lead  is  unlikely  to  be  isotopically  representative  of  the  whole  rock 
sample.  This  inhomogeneity  was  also  reported  by  Doe  and  Tatsumoto  ( 
1972)  for  lunar  samples.  Similar  results  were  also  observed  by  Sinha 
(1969)  and  Starik  et  al  (1957)  even  in  a  single  mineral  fraction  of 
the  rock.  Cumming  et  al  (1970)  cautioned  that  the  volatilization 
method  must  be  used  with  care  when  total  lead  composition  is  required. 
Black  et  al  (1971)  made  several  analyses  of  whole  rock  lead  composi¬ 
tion  both  by  volatilization  and  by  chemical  dissolution.  Among  the 
six  rocks  they  analysed,  three  rocks  were  reported  to  yield  consistent 


":>3  ■ 


isotopic  composition  by  both  methods.  However  in  one  rock  which  was 
said  to  be  more  radiogenic,  they  reported  a  gross  analytical  disagree¬ 
ment  of  about  6%  between  the  the  two  methods  with  the  volatilised 
lead  being  more  radiogenic  than  the  dissolution  lead.  They  suggested 
the  incomplete  lead  recovery  by  volatilization  as  a  reason  for  the 
discrepancy. 

As  mentioned  previously,  the  refractory  minerals  such  as  zircon 
contain  very  radiogenic  lead  and  therefore  the  isotopic  composition 
of  the  volatilised  lead  for  a  complete  recovery  should  be  more  radio¬ 
genic  than  that  of  slightly  incomplete  recovery.  It  is  difficult  to 
judge  whether  or  not  the  chemical  dissolution  method  gives  a  complete 
lead  recovery,  as  Black  et  al  reported  in  the  same  paper  that  the 
volatilization  method  yields  104%  of  lead  compared  to  the  chemical 
dissolution  method  reported  by  Tatsumoto  (1969)  which  was  made  on 
U.S.G.S.  standard  basalt  BCR-1. 

These  data  appear  to  indicate  that  the  volatilization  technique 
can  yield  complete  recovery  as  long  as  a  temperature  of  1300°C  or 
higher  is  attained. 

The  plots  of  Pb20^/Pb20^  versus  Pb20^/Pb20^  for  leads  volatilised 
at  different  temperatures  for  the  eight  sample  rocks  are  shown  in 
Figures  19  to  26.  Data  of  each  rock  fit  reasonably  well  on  a  single 
straight  line  (called  a  temperature  isochron  here)  as  expected  by 
Cooper  and  Green  (1969)  and  Cumming  et  al  (1970).  The  slopes  of  the 
best-fit  lines  through  the  data  points  were  calculated  by  the  method 
of  Cumming  et  al  (1972)  which  includes  the  effects  of  correlated 
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errors  in  the  two  variables.  For  most  of  the  analyses  the  sum  of  the 
squared  residuals  exceeds  the  number  of  degrees  of  freedom  indicating 
that  the  data  do  not  fit  the  line  within  the  experimental  error 
assigned  to  the  data.  The  discrepancy  may  indicate  that  there  is  some 
lack  of  closed  system  behaviour  for  the  rock  or  that  the  precision  of 
the  analyses  has  been  over  estimated.  As  a  first  approximation  the 
data  are  treated  as  if  they  fit  the  calculated  line  and  the  error 
estimates  of  the  slope  have  been  increased  by  the  ratio  (^residuals/ 
degrees  of  freedom)^.  This  procedure  has  the  effect  of  combining  the 
measuring  error  with  possible  geological  variation  but  at  the  present 
stage  of  this  work  it  is  believed  that  this  procedure  leads  to  the 
most  satisfactory  interpretation  of  the  data. 

The  linear  relationship  of  ratios  derives  from  the  fact  that  if 
the  rock  as  a  whole  contained  initial  amounts  of  lead  which  were 
homogeneous  in  composition,  then  the  lead  contained  in  different 
mineral  phases  would  changes  isotopically  with  time  according  to  their 
respective  U/Pb  ratios  and  at  any  particular  temperature  we  measure 
some  linear  combination  of  contributions  from  the  various  minerals. 
Thus  if  the  system  remains  closed,  the  slope  of  the  straight  line  may 
be  interpreted  in  terms  of  the  time  since  the  formation  of  the  system. 

6.3.  Two  stage  model  temperature  isochron  ages 

The  slope-ages  and  calculated  values  of  p  based  on  a  two-stage 
lead  evolution  model  are  shown  in  Table  13*  An  inspection  of  the  data 
shows  that  samples  YKWG  #3,  #4,  #7  and  #8  have  ages  of  about  2550  m.y. 
The  sample  YKWG  #6  which  is  a  pegmatite  dike  should,  according  to  the 
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Table  13 .  Slope  ages  of  the  temperature  isochrons 


Sample : 

p  value 
(95f«  C.L.) 

Slope 
(S.D. ) 

Age  in  m 
(S.D.) 

YKWG  #1 

8.60 

.113 

1880 

(Syenite,  acidic) 

±.02 

±.005 

±75 

YKWG  #2 

8.81 

.176 

2650 

(Syenite) 

+  .02 

±.036 

±300 

YKWG  #3 

8.65 

.168 

2575 

( Gr anodiorit e ) 

+  .02 

±.004 

±40 

YKWG  #4 

8.63 

.165 

2550 

(Gr anodiorit e) 

+  .04 

±.030 

±270 

YKWG  #5 

8.65 

.163 

2525 

(Gr anodiorit e) 

+  .02 

±.005 

±50 

YKWG  #6 

8.91 

.174 

2635 

(Pegmatite) 

±.03 

±.001 

+  15 

8.97 

.175 

2645 

±.02 

±.001 

.  ±15 

YKWG  #7 

8.65 

.158 

2474 

(Gr anodiorit e) 

±.01 

±.004 

+45 

YKWG  #8 

8.82 

.163 

2530 

(Diorite) 

±.02 

±.009 

±90 

geological  relationships,  be  younger  than  its  host  rock.  The  older  age 
suggested  by  the  lead  isotope  ratios,  may  for  the  time  being  be 
explained  on  the  basis  of  B  type  anomalous  lead  as  defined  by  Houtermans. 
Possibly  some  addition  of  radiogenic  lead  from  the  surrounding  rocks  to 
the  pegmatite  might  have  occurred  during  the  formation  of  the  dyke. 

The  age  of  YKWG  #1  (i860  m.y.  -  75  m.y.)  is  identical  within  experimen¬ 
tal  error  with  the  data  of  Gates  and  Hurley  (1973)  for  the  high  Rb 
samples  of  dyke  set  I.  In  fact  a  recalculation  of  the  age  of  these 
data,  using  5  samples  containing  the  highest  Rb/Sr  ratios,  yield  an 
age  of  1870  m.y.  +  50  m.y.  in  excellent  agreement  with  the  data  given 
here.  There  seems  to  be  no  reason  to  ignore  these  interesting  data  on 
the  diabase  dykes. 

6.4. _ Mineral  and  whole  rock  isochrons 

When  all  ratios  except  those  of  YKWG  #1  are  used  for  the  straight 
line  fitting  (Figure  27),  a  line,  called  "all-temperature"  isochron 
here,  with  slope  0.175  ±  .002  is  obtained  and  the  slope-age  is  2650  m.y. 
+  20  m.y.  intermediate  between  the  whole  rock  and  heating  results.  If 
the  temperature  isochrons  are  to  be  interpreted  as  the  equivalent  of 
mineral  isochrons  then  one  would  expect  this  relationship  and  no 
particular  significance  may  be  attached  to  the  slope  obtained. 

The  difference  in  ages  between  the  whole  rock  isochron  and  the 
average  of  the  individual  temperature  isochrons  may  be  significant  and 
the  time  interval  compares  well  with  the  period  of  tectonic  activity 
deduced  by  Green  and  Baadsgaard  (l97l). 

The  comparison  of  the  whole  rock  lead  isochron  with 
mineral  or  heating  temperature  isochrons  is  analogous  to  the 
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similar  situation  (Ulrych  and  Reynolds,  1966)  in  Rb/Sr  systematics 
which  are  well  known  (for  example,  see  Baadsgaard  1965)*  Let  the 
subscript  0  denote  the  time  at  which  a  rock  or  mineral  first  became 
a  closed  system  and  p  denote  the  present  time.  Combining  the  known 
relations 

(Rb87)0  =  (Rb87)/* 

and  (Rb87)Q  “  (RbS?)p  =(Sr87)p  (radiogenic  only) 

we  obtain  (Rb^7)  (exb  -  l)  =  (Sr^7)  . 

P  P 


Considering  the  change  in  the  total  amount 


of  Sr87, 


(Sr87)p  =(Sr87)o  -  (Sr87); 


we  obtain  the  equation 


Equation  (A)  is  the  equation  of  a  straight  line  as  shown  in  Figure  28 
for  samples  which  have  same  initial  (Sr  /Sr  )Q  ratios.  When  the  rock 
is  metamorphosed  at  time  tm,  the  Sr  in  the  rock  is  redistributed  and 
the  Rb/Sr  in  the  recrystallized  minerals  need  not  be  the  same  as 
formerly  so  the  slopes  of  the  growth  lines  may  be  different  after  meta¬ 
morphism.  The  whole  rock  growth  line  is  unaffected,  because  of  the  fact 

87  86 

that  the  rock  is  a  closed  system.  The  Sr  /Sr  in  the  recrystallized 

87  86 

(or  newly- formed)  minerals  grows  again  from  the  Sr  /Sr  values  for 
the  whole  rock  at  tm.  Figure  29  shows  the  whole-rock  and  mineral 
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Effect  of  metamorphism  on  the  change  m  Sr  /Sr 
with  time  in  a  rock  (closed  system)  and  its 
constituent  minerals  (open  systems). 

(after  Baadsgaard,  1965) 
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®  Whole  rocks 
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Figure  29  • 


A  —  whole  rock  composition  at  t  . 

B  —  whole  rock  composition  at  tm=R^  at  tm, 

the  new  initial  ratio  for  the  minerals 

at  the  time  t  . 

m 

C  —  whole  rock  composition  at  t  . 

D„  —  mineral  composition  at  t  after  resetting 
1-n  ^  p 

of  R  at  time  t  . 

m 

Variation  in  Sr^/Sr^  vs.  Rb^/Sr^  for  whole- 
rock  and  the  effect  on  constituent  minerals  by 
a  met amorphic  event  at  time,  tm» 


■ 


isochrons  for  two  whole  rock  phases  and  the  result  of  a  metamorphic 


event  at  t  . 

m 
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In  the  case  of  a  U-Pb  system,  the  equation  for  Pb  may  be 


written  as 

( 


p,  206  -206  tt238  x, 

Pb  L  -  -  (-^X*  *  -1) 


pb2<Vp  Vpb  20/X° 


Pb 


•(B); 


207  208 

equivalent  expressions  exist  for  Pb  and  Pb 

A  lead-lead  diagram  showing  the  whole  rock  isochron  and  the 
mineral  isochrons  similar  to  Figure  29  will  also  occur  for  lead 
isotopes  and  is  shown  in  Figure  30* 

The  difference  in  ages  between  the  whole  rock  isochron  and  the 
average  individual  temperature  isochrons  may  in  fact  represent  two 
different  ages.  A  whole  rock  isochron  age  (2700  m.y.  +  85  m.y.)  is 
the  time  when  the  rock  crystallized.  The  mineral  isochron  age 
2550  m.y.  ±  50  m.y.  is  equivalent  to  the  temperature  isochron  age, 
the  time  which  signals  the  end  of  a  metamorphic  event  when  the 
individual  mineral  phases  are  able  to  retain  their  radiogenic  lead. 


6. 5* _ Some  implications 

From  the  volatilized  lead  ratios  discussed  in  this  Chapter, 
some  important  implications  regarding  the  volatilization  technique 
and  its  application  to  geochronology  may  be  derived. 

(l).  The  temperature  isochron  obtained  from  the  incremental  heating 
of  a  rock  is  analogous  to  the  mineral  isochron  of  the  rock  as  claimed 
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Figure  30  .  Lead-lead  diagram  for  whole  rocks  and  minerals 
showing  a  two-event  history. 


by  Cumming  et  al  (1970). 

(2) .  The  all-temperature  isochron  age  obtained  from  the  incremental 
heating  of  a  number  of  cogenetic  rocks  is  intermediate  between  the 
whole  rock  isochron  and  the  mineral  isochron  ages. 

(3) .  Extraction  of  lead  from  a  number  of  cogenetic  rocks  by  incremen¬ 
tal  heating  of  each  rock  is  a  fast  and  efficient  method  of  finding 
some  information  about  the  history  of  the  rocks. 

To  test  the  first  suggestion  that  the  temperature  isochron  is 
analogous  to  the  mineral  isochron,  mineral  separations  were  made  on 
three  rock  samples  and  lead  from  each  of  the  mineral  phases  were 
extracted  by  volatilization.  Results  of  the  mineral  leads  are  discuss¬ 
ed  in  Chapter  VII. 

6.6.  Temperature  isochron  of  southeast  granodiorite 

As  mentioned  in  Chapter  I,  one  rock  sample  from  the  southeast 
granodiorite  obtained  from  the  Department  of  Geology  was  analysed. 
Table  14  shows  the  ratios  of  lead  extracted  at  various  different 
temperatures  and  the  lead-lead  diagram  is  shown  in  Figure  31  The 
temperature  isochron  has  a  slope  of  0.1767  ±  .0043(S.D.)  which  repre¬ 
sents  a  slope-age  of  2660  m.y.  ±  45  m.y..  Thus  if  this  sample  is 
representative  of  the  material  of  the  southeast  granodiorite,  it  then 
appears  that  the  mineral  age  of  the  southeast  granodiorite  is  similar 
to  that  of  the  western  granodiorite.  A  similar  relationship  was  also 
observed  by  Green  (1968)  using  the  Rb/Sr  decay  scheme. 
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Table  14.  Ratios  of  the  volatilized  lead  of  the  southeast 
granodiorite 


Volatilization 
t emper atur e ( 0 C ) 

Pb206/Pb204 

Pb207/Pb204 

Pb208/Pb20Z+ 

900 

21.738 

16.071 

43 . 177 

.020 

.015 

.040 

950 

20.701 

15.871 

41.568 

.062 

.048 

.126 

1000 

20.627 

15.9U 

41.619 

.023 

.018 

.047 

1050 

18.864 

15.608 

38.443 

.051 

.043 

.106 

1075 

18.584 

15.493 

38.133 

.024 

.020 

.049 

1100 

17.287 

15.303 

36.451 

.026 

.024 

.057 

17.240 

15.300 

36.381 

.066 

.059 

.140 

*  Internal 

estimate  of  standard 

deviation. 

TEMPERATURE  ISOCHRON 
SOUTHEAST  GRANODIORITE 
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CHAPTER  VII 


MINERAL  LEAD  COMPOSITION  AMD  THE  MINERAL  ISOCHRON 

7.1.  Introduction  and  mineral  separation 

To  test  the  prediction  made  by  Cumming  et  al  (1970)  that  the 
temperature  isochron  is  analogous  to  a  mineral  isochron,  mineral 
separations  were  made  on  three  samples;  YKWG  #1,  YKWG  #4  and  YKWG  #5* 
Mineral  separations  were  performed  with  standard  heavy  liquid 
techniques  involving  the  use  of  tetrabromo ethane,  methylene  iodide 
and  Clerici  solution.  The  analytical  procedures  are  described  below. 

(1) .  The  rock  was  first  broken  into  small  pieces  and  weathered  mate¬ 
rial  discarded.  The  remaining  material  was  then  crushed  to  about 
120  mesh. 

(2) .  A  Wilfley  table  was  then  used  to  separate  the  rock  powder  into 
three  gross  fractions,  the  light  (called  gravity  up  and  denoted  by 
W  f  ),  the  heavy  (called  gravity  down  and  denoted  by  W  I  )  and  the 
intermediate  (called  gravity  float  and  denoted  by  W-»  ). 

(3) .  After  drying  the  powder,  the  irons  (from  the  crusher)  were  picked 
out  with  a  hand  magnet.  Magnetic  separations  were  then  made  on  the 
light  fraction  (for  biotit e  and  feldspar)  and  the  heavy  fraction 

(for  apatite,  zircon,  pyrite,  ...  etc.).  Nonmagnetic  material  (denoted 
by  M-T  )  was  separated  from  the  magnetic  fraction  (denoted  by  Mi  ). 

(4) .  The  heavy  liquids  tetrabromoethane  (TBr^)  and  methylene  iodide 
(MI2)  were  then  used  for  density  separation. 

Because  of  the  low  content  of  zircon  and  the  small  sample  of 
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rock  available,  a  final  separation  of  zircon,  pyrite,  sphene,  ..  etc. 
was  not  possible  and  the  final  residual  mixture  was  treated  as  a 
mixed-mineral  phase  for  lead  extraction. 

The  combinations  for  the  separations  of  the  four  mineral  phases 
used  for  lead  extractions,  biotit e,  potassium  feldspar,  apatite  and 
the  mixture  of  zircon,  pyrite,  ...etc.  are: 

Biotite:  Wt  ,  Ml  ,  TBr^j 
Biotit e  impurity:  Wt  ,  Ml  ,  TBr^t 
K  feldspar:  Wt  ,  Mt  ,  dilute  TBr^t 
Apatite:  Wl  ,  Mt  ,  TBr^l  ,  MI^I 

Mixture  of  zircon,  pyrite,  etc.:  Wl  ,  Ml  ,  TBr^t  ,  MI^I 


7.2.  Mineral  lead  composition  and  the  mineral  isochron 

Table  15  shows  the  mineral  lead  composition  of  the  three  rock 
samples.  The  leads  were  extracted  by  induction  furnace  heating  and 
purified  by  the  dithizone  extraction  technique.  Ratios  shown  are 
normalized  to  N.B.S.  Equal  Atom  lead  ratios  given  by  Stacey  et  al 
(1969).  Each  of  the  minerals  was  heated  in  the  furnace  at  1350°C  for 
about  one  hour.  It  can  be  seen  that  lead  contained  in  K-feldspar  is 
least  radiogenic  followed  by  lead  in  biotite,  apatite  and  the  mixture 
of  zircon,  pyrite,  etc.. 

The  lead-lead  diagrams  are  shown  in  Figures  32,  33  and  34*  The 
mineral  lead  ratios  of  each  sample  fit  on  a  straight  line  with  a 
slope  which  agrees  within  the  standard  deviation  to  the  slope  of  its 
respective  temperature  isochron. 
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Table  15.  Mineral  lead  ratios 


Sample : 

Mineral  phase 

pb2°6/pb204 

Pb207/Pb20i 

Pb208/Pb204 

YKWG  #1 

K- feldspar 

14.975 

.010* 

15.207 

.010 

34.355 

.023 

14.930 

.014 

15.143 

.014 

34.259 

.034 

Biotite 

21.310 

.022 

15.907 

.016 

37.444 

.039 

21.298 

.036 

15.870 

.027 

37.331 

.066 

Biotite  impurity 

18.248 

.020 

15.517 

.020 

36.523 

.033 

Mixture  of  zircon, 
pyrite,  etc. 

no 

lead 

no 

lead 

no 

lead 

YKWG  #4 

Biotite 

_ 

— 

— 

K-feldspar 

14.169 

.019 

14.818 

.020 

35.004 

.052 

14.265 

.015 

14.905 

.016 

35.196 

.039 

Apatite 

87.981 

.145 

27.374 

.045 

90.899 

.153 

Mixture  of  zircon, 
pyrite,  etc. 

123.054 

.236 

34.159 

.066 

79.521 

.156 

122.835 

.501 

34.109 

.140 

79.442 

.328 

*  Internal  estimate  of  standard  deviation. 


-  D'« 

Table  15  cont'd 


Sample : 


YKWG  #5 


Mineral  phase 


K- feldspar 


Biotit e 
Apatite 


Mixture  of  zircon, 
pyrite,  etc. 


Pb206/Pb204 

Pb207/Pb2C* 

Pb2°8/Pb204 

15.608 

14.970 

35.621 

.013* 

.012 

.029 

15.582 

14.964 

35.596 

.013 

.013 

.031 

45.710 

19.019 

47.961 

.370 

.157 

.406 

45.482 

18.984 

47.645 

.375 

.150 

.381 

45.863 

19.088 

48.416 

.163 

.069 

.181 

135.467 

34.688 

63.990 

.279 

.072 

.134 

*  Internal  estimate  of  standard  deviation. 
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MINERAL  ISOCHRON 


114 


U-\ 

CV 


CP 

vO 

to 
\ — 1 

• 

rO 

to 

MO 

S 

-4 

CN2 

!>- 

o~\ 

• 

UA 

MO 

* 

O 

• 

• 

w 

O 

^ — 1 

-4 

to 

o 

II 

II 

II 

Ji 

o 

o 

o 

O 

PL 

Ph 

ctf 

rd 

-P 

0*9l 

*1QZ^  LOZ^ 


o 

cv 

C\2 


u~\ 

'O 

CNJ 


o 

'O 


-4 

o 

c\ 

-£2 

Ph 


'O 
lt\  O 
•  C\2 

■I>-  jO 

T— I  Ph 


o 

KvO 


•fr-4 


O 


- J- 

0*81 


- H 

9*91 


— h- 

9*  Cl 


0*21 


Figure  32 


115 


U"\ 


\ 


O 
-  0s 


tr\ 

O 

cv 


-4 

o 

C\2 

PL, 


4) 

&  £ 


4) 

T — I 


Lf\ 

■-4 


O 

cr\ 


0*81 


£*91  0*£T 

noz^Loz^ 


S’£t 


0*21 


Figure  33 


MINERAL  ISOCHRON 


116 


UA 


-3’ 

C<A 

<D 

U 

a 

id 


Table  16. 

Slope-ages  and p  values 

of  the  mineral 

isochrons 

Sample : 

Calculated  j i 
value  (95%  C.L. ) 

Slope  (S.D.) 

Age  in  m 
(S.D.) 

YKWG  #1 

8.66  ±  .02 

0.112  ±  .005 

1865  ±  85 

YKWG  #4 

8.67  ±  .03 

0.171  ±  .002 

2600  ±35 

YKWG  #5 

8.30  ±  .02 

0.162  +  .004 

2520  ±  40 

Slope  ages  and  p.  values  calculated  on  the  basis  of  a  two-stage 
lead  evolution  model  for  the  three  mineral  isochrons  are  listed  in 
Table  16.  The  parameters  a  ,  bQ,  and  tQ  used  for  such  calculations 
are  the  same  as  those  used  for  the  temperature  isochrons. 

The  results  of  mineral  and  temperature  isochrons  indicate  that 
both  isochrons  yield  a  consistent  slope-age  for  a  rock.  If  this 
consistency  is  general  (for  all  types  of  rocks),  the  lead  isotope 
dating  procedure  can  be  greatly  simplified  by  employing  the  volati¬ 
lization  technique  since  no  mineral  separation  is  required. 


7.3.  Sulfide  mineral  isochron  of  the  Yellowknife  volcanics 

To  study  the  relationship  in  ages  between  the  western  granodio- 
rite  and  the  Yellowknife  volcanics,  twenty-three  sulfide  and  five 
galena  samples  separated  mostly  from  the  Yellowknife  volcanics  (two 
sulfides  and  one  galena  from  the  Yellowknife  sediments)  were  analysed. 
Sample  collection,  mineral  separations  and  lead  extractions  on  these 
samples  were  made  by  A.  P.  Leech  of  the  Department  of  Geology.  The 
writer  made  the  mass  spectrometric  analyses  and  the  results  are  shown 
in  Tables  17  and  18. 

Figure  35  shows  the  plot  of  ratios  Pb20^/Pb20^  versus  Pb20^ /Pb20^ 
of  the  twenty-three  sulfide  samples.  The  data  show  very  much  scatter 
around  the  best  fitting  straight  line  which  has  a  slope  of  0.173  ± 
.0064.  The  scatter  is  believed  mainly  due  to  various  sources  of 
analytical  error.  These  analyses  were  performed  during  the  initial 
stage  of  this  work  and  the  precision  is  not  as  good  as  that  obtained 
later.  The  ji  value  is  calculated  to  be  8. 90  ±  .04. 
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Table  17.  Pb  ratios  of  sulfide  samples  from  the 
Yellowknife  volcanics 


Sample  ident . : 

pb206/pb20^ 

pb207/pb204 

Pb2°8/pb204 

L-9,  pyrite 

13.843 

.013* 

14.801 

.014 

33.386 

.032 

L-8,  pyrite 

13.870 

.020 

14.858 

.021 

33.522 

.048 

L-8(+G) ,  chalcopyrite 

13.892 

.058 

14.871 

.062 

33.580 

.141 

L-7,  pyrite 

13.900 

.009 

14.885 

.010 

33.619 

.022 

L-8(-G),  pyrite 

13.901 

.010 

14.893 

.011 

33.622 

.026 

13.909 

14.898 

33.644 

Ptarmigan  pyrite 

14.017 

.012 

14.947 

.013 

33.675 

.026 

14.033 

.014 

14.964 

.015 

33.718 

.035 

L-26,  pyrite 

14.131 

.041 

14.971 

.044 

33.795 

.100 

L-29,  arsenopyrite 

14. 164 
.036 

14.817 

.037 

33.541 

.084 

14.173 

.061 

14.879 

.064 

33.646 

.  146 

Giant  15-8,  pyrite 

14.136 

.066 

15.058 

.071 

34.096 
.  166 

L-19,  pyrite 

14.282 

.013 

15.026 

.013 

34.347 

.032 

Ptarmigan  chalcopyrite 

14.596 

.051 

15.072 

.052 

34.675 

.121 

*  Internal  estimate  of  standard  deviation. 
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Table  17  cont'd 


Sample  ident . : 

Pb206/Pb204 

pb2°7/pb204 

pb208/pb204 

K-97-117,  pyrite 

14.734 

.029* 

15.110 

.030 

34.790 

.069 

L-6,  pyrite 

14.867 

.034 

15.188 

.035 

34.554 

.081 

K-97-116,  pyrite 

14.990 

.045 

15.150 

.046 

34.797 

.106 

L-21,  pyrite 

15.312 

.041 

15.317 

.042 

35.697 

.102 

B-1143,  pyrite 

15.814 

.023 

15.225 

.022 

35.508 

.051 

L-2,  pyrite 

15.977 

.030 

15.532 

.029 

34.430 

.067 

15.980 

.036 

15.526 

.035 

34.392 

.082 

L-24,  arsenopyrite 

16.341 

.089 

15.475 

.084 

36. 424 
.204 

Crestaurum  chalcopyrite 

17.015 

.026 

15.414 

.023 

36.131 

.055 

B-1269,  pyrite 

17.476 

.020 

15.468 

.018 

37.747 

.044 

Akaithcho  chalcopyrite 

18.277 

.016 

15.687 

.014 

38.235 

.035 

Akaithcho  pyrite 

17.937 

.011 

15.530 

.010 

37.626 
.  024 

L-20,  pyrite 

22.094 

.042 

l6 . 426 
.031 

49.947 

.095 

22.119 

.063 

16.441 

.047 

49.796 

.143 

*  Internal  estimate  of  standard  deviation. 


Table  18.  Pb  ratios  of  galena  samples  from  the  Yellowknife 


volcanics 


Sample  ident.: 


Pb206/pb2Q/j-  Pb207/Pb2°4  Pb2Q8/Pb2Q/f 


Ptarmigan  galena 

13  .844 
.030* 

14.891 

.032 

33.773 

.075 

13.870 

.028 

14.915 

.031 

33.846 

.070 

L-8  Rycon  galena 

13.871 

.012 

14.852 

.013 

33.536 

.031 

Con  galena,  3300 
level-Cambell  shear 

13.934 

.019 

14.978 

.021 

33.820 

.048 

Con  galena 

13.937 

.016 

14.977 

.017 

33.818 

.039 

Rycon  galena,  3300 
level-Campbell  shear 

13.944 

.014 

14.971 

.015 

33.870 

.034 

*  Internal  estimate  of  standard  deviation. 
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The  spread  of  ratios  along  a  straight  line  is  believed  to  be 
the  result  of  less  than  perfect  mineral  separations,  and  the  lead 
ratios  observed  are  not  representative  of  leads  in  pyrites  (assorted 
sulphides)  but  are  linear  combinations  of  contributions  from  sulfides 
and  different  proportions  of  impurities  remaining  after  mineral 
separation  (possibly  zircon,  sphene,  etc.)*  Microscopic  examination 
of  some  of  the  samples  indicated  that  various  impurities  were  present 
(Cumming  personal  communication,  1973)  and  in  many  cases  separation 
was  extremely  difficult  due  to  the  great  range  of  minerals  present# 

The  best  fitting  straight  line  thus  represents  a  mineral  isochron  of 
the  volcanic s  or  at  least  the  mineralized  rock  in  the  ore  zones. 

Although  the  observed  ratios  are  not  representative  of  the 
sulfides,  an  age  can  still  be  derived  from  the  mineral  isochron  and  is 
found  to  be  2630  m.y.  ±  70  m.y.  from  the  slope  of  the  line  given 
above.  If  we  take  the  least  radiogenic  lead  from  the  sample  L-9  to 
represent  the  lead  in  sulfide,  the  ratio  lies  on  the  conformable 
growth  curve  at  2630  m.y.  identical  to  the  slope— age  of  the  mineral 

isochron. 

Figure  36  shows  the  lead-lead  diagram  of  the  five  galena  samples. 
The  ratios,  which  agree  with  each  other  within  0 .5#  lie  slightly  above 
the  conformable  growth  curve  at  26i+0  m.y.,  the  same  age  as  the  pyrite. 
Thus  the  mineralization  of  galenas  and  pyrites  (metallic  sulphides)  of 
the  Yellowknife  volcanics  occurred  contemporaneously  at  about  263O  m.y. 
+  70  m.y.  ago.  The  volcanics  should  then  be  at  least  this  old. 
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7,4.  An  indication  of  the  age  of  the  Yellowknife  metasediments 

Among  the  twenty— three  sulfide  and  five  galena  samples  analysed 
two  sulfides  (Ptarmigan  pyrite  and  Ptarmigan  chalcopyrite)  and  one 
galena  (Ptarmigan  galena)  are  separated  from  the  Yellowknife  meta¬ 
sediments.  Lead  ratios  from  the  two  sulfide  samples  are  not  interpre 
table  due  to  possible  impurities  from  mineral  separations.  The 
Ptarmigan  galena  has  ratios  equal  to  those  of  the  galenas  from  the 
volcanics  which  suggests  that  the  age  of  the  sediments  should  at 
least  be  as  old  as  the  galenas  at  about  2640  m.y.. 
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CHAPTER  VIII 


SUMMARY  AND  CONCLUSIONS 


8.1.  Introduction 

It  has  been  pointed  out  by  many  workers  that  the  U/Pb  system 
has  distinct  advantages  over  other  methods  of  dating  because  one  can 
obtain  useful  information  from  the  Pb  isotope  ratios  alone  without 
the  necessity  of  chemical  extraction  of  parent  and  daughter  elements. 
In  this  work  the  advantages  of  the  Pb/Pb  method  have  been  extended 
so  that  it  is  possible  to  obtain  both  whole  rock  and  mineral  isochron 
information  on  the  same  rock  in  a  simple  manner.  The  extra  informa¬ 
tion  can  be  interpreted  geologically  in  a  manner  consistent  with 
age  measurements  by  other  techniques, 

8.2.  Interpretation  and  comparison  with  other  results 

A  whole  rock  lead-lead  isochron  age  of  2700  m.y.  ±  85  m.y.  was 
derived  from  six  western  granodiorite  samples  based  on  a  two  stage 
model.  Since  the  concordia  plot  of  Figure  17  indicates  that  a  two 
stage  model  cannot  account  for  the  observed  ratios,  it  is  necessary 
to  consider  those  three  stage  models  which  yield  rectilinear  plots. 
For  our  data,  the  best  straight  line  passes  through  the  intersection 
of  the  zero  isochron  and  the  conformable  growth  curve,  so  if  we 
interpret  the  data  in  terms  of  a  three  stage  f ^constant  model,  the 
lower  intercept  must  correspond  to  t2  and  therefore  there  would  have 
been  an  earlier  event  at  the  time  t ^ ,  with  an  age  greater  than  2700 
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million  years.  There  is  no  support  in  the  present  data  for  such  an 
event  although  Robertson  and  Cumming  (1968)  obtained  a  lead-lead 
model  age  of  about  4000  m.y.  based  on  analyses  of  some  anomalous 
leads  from  galena  samples. 

The  other  three  stage  model  which  yields  a  linear  lead-lead  plot 
is  the  case  where  f^constant  for  which  the  lower  intercept  yields  t^ 
regardless  the  value  of  f^  and  it  is  impossible  to  test  this  parti¬ 
cular  model  any  further  from  the  lead-lead  data  alone. 

From  the  concordia  plot  of  Figure  17,  the  best  line  through  the 
data  points  intersects  the  normal  concordia  curve  at  1540  ±  200  m.y. 
and  if  a  three  stage  model  applies,  this  time  is  an  upper  limit  to 
t^  for  the  f  ^constant  case  or  equal  to  t^  if  f^=  const  ant,  i.e. 
t  <  1540  ±  200  m.y..  The  only  indication  of  an  event  at  this  time 
is  a  poorly  defined  Rb/Sr  isochron  on  dyke  set  III.  Gates  and  Hurley 
(1973)  obtained  an  age  of  1370  m.m.<t<  1730  m.y.  for  these  rocks. 

It  is  possible  that  the  intrusion  of  these  dykes  may  have  been  related 
to  the  end  of  a  second  stage  of  the  lead-uranium  system  in  the  western 
granodiorite. 

The  lead-lead  mineral  isochron  age  of  2530  m.y.  ±  20  m.y.  (S.D. 
of  the  weighted  mean  of  the  rocks  YKWG  #2,  #3,  #4,  #5,  #7,  and  #8) 
cannot  correspond  to  in  these  three  stage  models  since  t2<  1540  m.y. 
+  200  m.y.,  therefore  it  must  represent  some  event  which  is  not 
reflected  in  the  whole  rock  data. 

Figure  37  summarizes  the  pertinent  age  measurements  of  the 
major  geological  units  of  the  Yellowknife  area.  Upon  comparison  of  the 
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Figure  37.  Age  measurements  of  the  Yellowknife  area 
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data  shown  in  Figure  37,  some  of  the  features  noted  are; 

(1) .  The  lead— lead  whole  rock  isochron  age  of  the  western  granodio- 
rite  agrees  reasonably  well  with  the  Rb/Sr  age  obtained  by  Green  et 
al  (1968). 

(2) .  The  mineral  age  of  the  western  granodiorite  is  about  150  m.y. 
less  than  the  whole  rock  age.  This  difference  of  ages  agrees  with 
the  cooling  interval  of  140  m.y.  deduced  from  K/Ar  mica  ages  by  Green 
and  Baadsgaard  (1971)* 

(3) .  The  age  of  the  pegmatite  dyke  based  on  the  mineral  isochron 
(temperature  isochron)  is  less  than  the  whole  rock  age  of  the  country 
rock.  The  age  of  the  pegmatite  dyke  apparently  has  not  been  updated 
by  the  process  which  altered  the  country  rock  at  2530  m.y.. 

(4) .  Lead  data  from  fractional  heating  of  a  southeast  granodiorite 
rock  sample  confirmed  the  Rb/Sr  analyses  that  this  granodiorite  is 
older  than  the  western  granodiorite. 

(5) .  The  much  younger  mineral  age  of  YKWG  #1  agrees  remarkably  well 
with  the  Rb/Sr  data  from  Gates  and  Hurley  (1973)  on  the  high  Rb  rocks 
of  dyke  set  I  (age  and  error  shown  in  Figure  37  are  recalculated 
values  using  five  rocks  which  have  the  highest  Rb/Sr  ratios). 

(6) .  The  mineralization  age  of  2630  m.y.  derived  from  galena  and 
and  pyrite  (metallic  sulphide)  samples  from  the  Yellowknife  group 
suggests  that  the  formation  of  volcanics,  intrusion  of  granite  and 
ore  formation  must  have  taken  place  in  a  short  period  of  time, no  more 
than  a  few  tens  of  million  of  years. 

(7)  .  An  apparent  conflict  in  age  relationships  appears  in  the  data 
of  Gates  and  Hurley  (1973)  for  "the  dyke  set  I.  The  dyke  which  they 
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sampled  cuts  the  Prosperous  Lake  granite  and  yet  appears  to  have  an 
older  age.  It  amy  be  that  the  Rb/Sr  mineral  isochron  age  on  the  granite 
corresponds  to  the  younger  lead-lead  mineral  age  and  that  the  rock  is 
actually  older. 

8.3 . _ Mass  discrimination 

Discrimination  in  a  mass  spectrometer  is  always  a  serious  problem 
in  precise  measurements  of  isotopic  abundances.  It  has  been  shown  that 
different  filament  positions  between  runs  in  a  single  filament  solid- 
source  mass  spectrometer  may  introduce  mass  dependent  errors  in  the 
isotope  ratios.  Some  evidence  of  possible  focusing  effects  on  the 
isotope  ratios  were  also  demonstrated. 

The  solution  to  this  problem  seems  to  be  to  ensure  that  the  source 
geometry  is  maintained  as  constant  as  possible  so  that  significant 
changes  of  focusing  parameters  are  not  required  between  runs. 
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APPENDIX  I 


SOME  FORMULAS  ON  VECTOR  DIFFERENTIATION 


T  T 

Let  x  and  y  be  column  vectors  and  x  ,  y  be  the  transpose  of 

T 

x  and  y  respectively,  let  A  be  a  square  matrix  and  A  be  the  transpose 
of  A.  Then, 


^(xTAy)=Ay 

9  /  T .  \  T.T 
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For  a  symmetric  matrix,  B,  we  have 

— |-(xTBx)=2Bx 
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152  CALL  TSONG (R, NVAR) 

153  CALL  MINV  (B,NVAR,FR,  L,M) 

154  IF (FR. EQ.  0, 0D0)  GO  TO  8 

155  GO  TO  9 

156  C  RHS  =  RIGHT  HAND  SIDE  VECTOR 
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182  NPl=NPEAK5-1 

183  DO  11  1=1, NP1 

184  11  R(I+1)=1.0D0/X(I) 
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APPENDIX  III 


DITHIZONE  EXTRACTION  PROCEDUES 


All  chemical  ware  was  initially  cleaned  with  distilled  water  and 
placed  in  hot  HNO^  for  several  hours  or  overnight  before  a  final 
rinsing  with  pure  water  before  being  used*  The  following  is  an  outline 
of  the  procedures  used  in  purifing  the  leads  for  mass  spectrometer 
analysis. 

(1) .  Adjust  the  pH  of  the  lead  containing  solution  of  3N  HNO^  with 
ammonium  hydroxide  (2 8%)  to  about  10  in  a  separatory  funnel. 

(2) .  Add  lc.c.  of  KCN  (in  2 %  NH^OH)  and  1  c.c.  of  30%  ammonia  citrate 
to  the  solution  to  form  complexes  of  other  reacting  metals. 

(3) .  Shake  the  solution  with  some  dithizone  ( 0.1 %  in  Chloroform).  The 
reaction  of  dithizone  with  lead  gives  a  brilliantly  colored  pink 
solution.  Transfer  the  chloroform-dithizone  solution  containing  lead 
into  a  smaller  separatory  funnel. 

(4) .  Repeat  step  (3)  until  the  solution  no  longer  shows  pink  color 
indicating  all  the  lead  has  been  taken  up. 

(5) .  Add  1.5  c.c.  of  2%  HNO^  into  the  pink  solution  and  shake  it,  the 
lead  is  redissolved  in  the  2%  HNCL.  Discard  the  dithizone. 

(6) .  Wash  the  2%  HNO^  solution  containing  the  lead  with  about  5  c.c. 

of  chloroform  at  least  twice  by  adding  the  chloroform  into  the  solution 
and  vigorously  shake  the  combined  solution.  Discard  the  chloroform. 

(7) .  Transfer  the  washed  2 Jo  solution  from  the  top  of  the  separatory 
funnel  into  a  2  c.c.  centrifuge  tube  and  adjust  the  pH  of  this  solution 
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to  about  4.5  and  bubble  it  with  H2S  gas  causing  PbS  sample  to  be 
precipitated. 

(8).  Centrifuge  the  precipitated  PbS  ready  for  loading  onto  the 
filament  for  mass  spectrometer  analysis. 


I 


B30090 


